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Abstract. We revisit the problem of Q-colourings of the triangular lattice using a mapping 
onto an integrable spin-one model, which can be solved exactly using Bethe Ansatz 
techniques. In particular we focus on the low-energy excitations above the eigenlevel g 2 , 
which was shown by Baxter to dominate the transfer matrix spectrum in the Fortuin- 
Kasteleyn (chromatic polynomial) representation for Qo < Q < 4, where Qo = 3.819 671 • • •. 
We argue that g 2 and its scaling levels dehne a conformally invariant theory, the so-called 
regime IV, which provides the actual description of the (analytically continued) colouring 
problem within a much wider range, namely Q € (2,4]. The corresponding conformal 
field theory is identihed and the exact critical exponents are derived. We discuss their 
implications for the phase diagram of the antiferromagnetic triangular-lattice Potts model 
at non-zero temperature. Finally, we relate our results to recent observations in the field of 
spin-one anyonic chains. 


1. Introduction 

The four-colour problem is probably the best known problem of graph theory, and has 
fascinated generations of mathematicians and laymen since it was hrst stated in a letter 
from Augustus de Morgan to Sir William Rowan Hamilton dated 23 October 1852. The 
suspense was resolved in late July 1976 when Appel and Haken announced the proof of 
the four-colour theorem: Every planar graph admits a vertex 4-colouring. 

A quantitative version of the problem was proposed by Birkhoff in 1912: Given a 
simple unoriented graph G and a set of Q colours, how many proper vertex colourings does 
G admit? This dehnes the chromatic polynomial Xg(Q)- The four-colour theorem can then 
be stated: If G is planar, then Xg(4) > 0. 
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Although initially defined for Q G N, it is easily proved j2]-[^ that Xg{Q) is in fact a 
polynomial in Q, and can as such be evaluated for any Q G C. In particular, if G = (V, E) is 
a simple unoriented graph with vertex set V and edge set E, its chromatic polynomial can 
be written as 

x'gW) = (1) 

ACE 

where the sum is over all spanning subgraphs (V, A) of G, and k{A) is the number of 
connected components of (V,A). This suggests an algebraic or even analytic approach to 
the colouring problem. There exist many studies of the location in C of the roots of Xg{Q), 
henceforth called chromatic roots. These studies concern either specific graphs, or all planar 
graphs, or other infinite families of graphs. 

The chromatic polynomial is a specialisation of the more general Tutte polynomial [^. 
Moreover, in statistical mechanics, the chromatic polynomial arises as a special case of the 


Q-state Potts-model partition function [m. (See section 2.1 


Many studies have aimed at elucidating the phase diagram of the chromatic polynomial, 
the nature of the different phases, and the transitions between them. Although a few results 
are available for non-planar graphs , most of the effort has concentrated on two dimensions 
{G planar or embedded in a cylinder or a torus), a case to which we specialise henceforth. 

Two lines of investigation have proven particularly powerful. In the first, the phase 
diagram is inferred from the study of chromatic roots for Q G C (see e.g. the series of papers 
[8|-[I^) for recursive graph families, corresponding to (finite or infinite) pieces of the square 
and triangular lattices with various boundary conditions. Similar but quantitatively more 
precise results have recently been obtained via a topologically weighted graph polynomial 


13 16 on more general families of graphs. 


In the second series of studies, the colouring problem has been exactly solved in the 
thermodynamic limit along the chromatic line Q G M for the triangular lattice [^[^. In 
this context, solvability of this model means that the exact free energy density (per site) 


/g(Q) = ^ logXG(Q) 


( 2 ) 


can be computed in closed form. On this solvable line, Xg{Q) can be transformed 
into an equivalent vertex model with a Yang-Baxter integrable i?-matrix, leading to an 
infinite family of commuting transfer matrices. Therefore, the triangular-lattice chromatic 
polynomial is integrable for any Q G M. 

The chromatic polynomial on the triangular lattice has only been studied directly 
through the computation of the bulk free energy 17, 18 , or indirectly via a series of 


transformations to an 0(n) loop model on the hexagonal lattice 19 . Baxter found 17 that 


in the complex Q-plane the free energy density ([^ is described by three distinct analytic 
functions gi{Q) {i = 1,2,3). In addition, he computed the phase diagram of the chromatic 
polynomial in the complex Q-plane (see figure [^, by finding on which region each function gi 
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Tri Lattice Chromatic Poly. 



Re(Q) 


Figure 1: Phase diagram of the triangular-lattice chromatic polynomial in the complex Q- 
plane . We show the three phases of the model; in each of them, one of the three functions 


gi{Q) found by Baxter 17,18 dominates, and the curves correspond to the loci where two 
of them coincide in absolute value: \gi\ = \gj\ > \gk\- At the bifurcation points A and B 
{Q ~ 3.859 627 ± 0.203 154z), all three functions are equimodular. The outer phase (where 
gi dominates) is non-critical. The other two phases are critical: the ground state of regime I 
(resp. regime IV) is given by g^ (resp. g 2 ). But one can extend analytically the ground 
state and its excitations of either of these two regimes inside the other, so that in each 
regime, we have the superposition of two distinct critical spectra. The colouring problem is 
actually described by the (analytically continued) ground state and excitations of regime IV. 
Furthermore, in regime I at Q = Sfc (|3), both the ground state g^ and its scaling levels will 
disappear from the partition function (^. 


was dominant, and locating the loci where two of these functions are dominant and coincide 
in absolute value 


18 


In this general picture, it is clear that the chromatic polynomial is critical only in the 
interval Q G [0,4]. Outside this interval, the model is non-critical, meaning that correlations 
of local observables decay exponentially fast with the distance. However, the determination of 
the critical exponents (which are expected to vary continuously, as in other integrable systems 
with a free deformation parameter—here Q), and the identihcation of the corresponding 
conformal held theory (CFT) are two important ingredients that are so far still missing. 

A major difference with other well-known integrable models is that the critical interval, 
Q G [0,4], allows for not just one, but two analytically unrelated expressions g 2 and g^ EZl 
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for the free energy ([^. Each one dehnes a different “phase” or regime. We shall call regime I 
the critical theory described by the ground state g^ and the set of conformally invariant 
scaling levels above it. Similarly regime IV will denote the critical theory corresponding to 
the ground state 5 ^ 2 • 

It is important to understand that this dehnition of regimes I and IV does not rely on any 
considerations whether the corresponding ground state eigenvalue—that is, for regime I, 
and g 2 for regime IV—is dominant in absolute terms. Baxter has indeed argued that 
g^ is dominant (i.e., the largest in absolute value among all transfer matrix eigenvalues) on 
the interval 0 < Q < Qo, with Qo = 3.819 671 ■ • •, whereas g 2 is dominant on Qo < Q < 4— 
but this is not the point here. We shall instead show that regime I makes sense for all 
0 < Q < 4, whereas regime IV makes sense (at least) for 2 < Q < 4. The scaling levels 
of either regime can in principle be followed (by analytic continuation) in the Bethe Ansatz 
(BA) computations, or by direct diagonalisation of the transfer matrix (TM), even into a 
range where the regime’s ground state is sub-dominant with respect to that of the other 
regime. So one of our principal hndings is that regimes I and IV are two superimposed 
critical theories that coexist, so to speak, as a pair of interchangable Russian dolls. 

This complication is presumably concomitant with the fact that the standard 
programme in statistical mechanics to hnd phase transitions and determine their critical 
properties in a given system, has not yet been carried out for this model. One of the main 
results of the present work is to identify the critical exponents and GET of both regimes. 
Interestingly, the more exotic regime IV is not present in the closely related hexagonal-lattice 
0(n) loop model, and it turns out to be related with a GET describing a spin-one anyonic 
chain (^ . 

Using different techniques, including both BA and TM computations, we moreover 
conclude that regime IV is the physical meaningful continuation of the critical colouring 
problem to Q G M, whereas regime I can best be described (from this point of view) as some 
universal and unphysical “junk” that has no relation to Q-colourings of the triangular lattice. 
In particular, we shall see that regime I—namely its ground state gs and the corresponding 
excitations—will disappear from the partition function ([^ at values of Q that are relevant 
to the colouring problem; i.e., the Beraha numbers 

Bk = 4 cos^ j , /c = 2,3,4,... . (3) 

In some sense, that will be made precise below, regime IV at the Beraha numbers therefore 
provides the maximal physical extension of the colouring problem with local observables. 
This is why a large part of this work is dedicated to unravelling the critical properties of 
regime IV. 

Before going on, we outline the contents of the paper. In section we will introduce 
the more general Q-state Potts model, which, in addition to the number of states or colours 
Q, has a temperature-like parameter v. In this section, we will describe the main features 
of the phase diagram in the real (Q,n)-plane of the Q-state Potts model on a 2D regular 
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lattice, using the best-known square lattice case as a guide. We will make emphasis on 
its antiferromagnetic (AF) and unphysical regimes. This section will provide the necessary 
background to discuss in detail the lesser known phase diagram of the triangular-lattice Potts 
model. 

After this introductory section, we will start to discuss our results. Section [^constitutes 
the hrst step in our programme. It provides a series of mappings that connect the triangular- 
lattice chromatic polynomial Xg{Q) with a spin-one vertex model on the square lattice. The 
former model is dehned on a strip graph with cylindrical boundary conditions. A key point 
in all these mappings is the introduction of a twist that ensures that the original boundary 
conditions are treated correctly. This ingredient was not essential to obtain the bulk free 
energies 17,^ (and hence was not discussed in those papers); but it has proven to be 


important in order to obtain critical exponents. 

Our Bethe-Ansatz work on the above-mentioned spin-one vertex model is presented in 
section [^ Results for regime I can be taken over, without much ado, from existing work on 
an integrable 7-vertex model on the honeycomb lattice, upon imposing the correct twist 
that ensures its equivalence with the triangular-lattice Potts model. The main endeavour 
however concerns regime IV, where two Fermi seas are present. In section [^ we obtain (using 
analytic techniques and some numerical checks) the conformal weights of various operators, 
culminating in the precise identihcation of the corresponding coset CFT for regime IV. 
The corresponding rational theories take the form of RSOS models, which we identify with 
recently studied spin-one anyonic chains 20 in section [^ Finally, section [^ contains our 
conclusions. 


2. Basic setup 

In this section, we will hrst describe the Q-state Potts model in the spin and Fortuin- 
Kasteleyn (FK) representations. In the AF case, the zero-temperature limit of the partition 
function specialises to the already dehned chromatic polynomial ([^. We will then discuss 
the physics of the best-known example: the square-lattice Potts model. Many features of 
the corresponding phase diagram are generic to all two-dimensional lattices, including the 
existence of an unusual phase called the Berker-Kadanoh phase. A more detailed discussion 
of this phase is contained in a separate subsection. Finally, we will review what is known 
about the triangular-lattice Potts model. 

2.1. The Q-state Potts model 

In statistical mechanics, the chromatic polynomial arises as a special case of the Potts 
model [^. Let G = {V, E) be an unoriented graph with vertex set V and edge set E, 
and attach to each vertex i E V a. spin variable cXi G {1,2,...,Q}, where Q G N. The 
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= E IT 

Wi} (*j> 


K&a a 


( 4 ) 


where K = J/T, the interaction energy between adjacent spins is — and T denotes 

the temperature. It is easy to see that in the AF (J < 0) case, the colouring problem arises 
in the T —0 limit, namely 


Xg{Q) = 


lim 


K^—oo 


Fortuin and Kasteleyn 
rewritten as: 

^^Potts 


22 


proved that the partition function 


Potts 


{Q,K) = = 

ACE 


yFK 


{Q,v) 


( 5 ) 

(|^ could be 

( 6 ) 


where |A| is the number of edges in the subset AGE, k{A) denotes the number of connected 
components in the spanning subgraph {V,A), and the temperature-like variable v is given 
by 

V = -1. (7) 

Therefore, the Potts-model partition function is a polynomial in the variables g,n, and the 
polynomial character of its chromatic specialisation, Xg{Q) = —1), is now obvious 

[cf. ([^]. This FK representation of the Potts model is useful to make sense of non-integer 
values of Q and/or imaginary values of the coupling constant J (i.e., v < —1). In this 
representation the AF regime corresponds to n G [—1,0), while the ferromagnetic one is 
given by n G (0, oo). 

As for the chromatic polynomial, many studies have aimed at elucidating the phase 
diagram of the Q-state Potts model, mostly in two dimensions. In particular, the zeros of the 


partition function have been investigated (see 23,24 and references therein) for recursive 
strip graphs of the square and triangular lattices with various boundary conditions. In 
addition, the use of the topologically weighted graph polynomial dehned in 13 16 , has 
provided more precise results on more general families of graphs. 

Remarks. 1. The Potts-model partition function (|^ can be dehned on any unoriented 
graph G, not necessarily simple. 

2. In graph theory the two-variable polynomial [cf. (|^] is known—after an 

innocuous change of variables—as the Tutte polynomial [sl. 


3. The authors of Refs. 23,24 considered the partition-function zeros for integer values 


of g > 2 (resp. real values of u > — 1) in the complex v (resp. Q) plane. 


2.2. The sguare lattice 

The phase diagram of the Potts model is best understood for the square-lattice model. In 


particular, Baxter 25,26 was able to solve exactly this Potts model in the thermodynamic 
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limit on specific curves in the (Q,n)-plane. On these solvable curves the model becomes 
integrable in the sense discussed in the Introduction. In particular, the square-lattice Potts 
model is integrable along the selfdual curves [25] 

( 8 ) 


= ^^/Q 


and along two mutually dual AF curves 26 


v± = 2 ± \/A — Q 


(9) 


Integrability studies lead hrst to the exact expression of the bulk free energy, from 
whose singularities some features about the phase transitions can be inferred. It is found 
that when moving across the self-dual curves ([^ in the u-direction, the system undergoes 
a hrst-order (resp. second-order) phase transition when Q > 4 (resp. when 0 < Q < 4). 
When crossing (§ the transition is simultaneously hrst and second order [^. Moreover, 


when moving across Q = 4 in the direction along the integrable curves, the free energy 
exhibits an essential singularity There are also singularities in the surface and corner 


free energies in this case 29,30 


Within the critical regime 0 < Q < 4 much further information can be obtained. Along 
the integrable curves one can compute various critical exponents, which will typically vary 
continuously with Q. These exponents provide crucial information that can help identifying 
the CFT describing the continuum limit. This program has been carried out for the self dual 
curves (|^, and the corresponding CFT is found to be that of a compactihed boson, that 
can be effectively described by the Coulomb gas (CG) technique 31 . This is arguably the 
simplest possible CFT with continuously varying exponents. 

The coupling constant g of the CG that describes the continuum limit along the selfdual 
curves (1^ is related to the temperature-like parameter v as: 


V = —2cos(7r5f), Q = 


( 10 ) 


This CG coupling constant satishes | < 1 for the ferromagnetic branch v = +\/Qi while 

0 < 5 ^ < I for the other (‘unphysical selfdual’ 32 ) branch v = —\/Q- The thermal operator 


has critical exponent 31 33 


XT 


^9 ^ 


( 11 ) 


and is conjugate to a perturbation in the temperature variable v around the critical curve 
(|^. In particular, the temperature perturbation is relevant {xt < 2) along the upper branch 
of ([^, which is identihed with the ferromagnetic critical curve vf{Q) = +y/Q- On the other 
hand, this temperature perturbation is irrelevant {xt > 2) along the lower branch of ([^, 
which will hence act as a basin of attraction for a hnite range of u-values. 

This basin is delimited by the AF curves ([^ and is called the Berker-Kadanoff (BK) 
phase 32 . The BK phase is however unphysical, in the sense that all the scaling levels 


corresponding to the GG description disappear from the spectrum whenever Q is equal 
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to a Beraha number Bf^ (|^. Two distinct representation theoretical mechanisms are 
responsible for this phenomenon. First, the multiplicity (also known as amplitude, or 
quantum dimension) of certain eigenvalues of the corresponding TM vanishes at Q = as 
can be seen from a combinatorial decomposition of the Markov trace Second, the 

representation theory of the quantum group t/q(s/ 2 ) for q a root of unity (with ^/Q = q + q“^) 


guarantees that other eigenvalues are equal in norm at Q = B^ 32,36 , and as their combined 
multiplicity is zero, they vanish from the spectrum as well. As a result, for Q = Bk only 


local observables remain, and these can be realised in an equivalent RSOS height model 37 


Finally, the continuum limit along the curves v± (§ that bound the BK phase, gives 
rise to a more exotic CFT 1^1^ with one compact and one non-compact boson that 


couple to form the SL(2,M)/U(1) black hole sigma model 39 , familiar in the string-theory 


context 40,41 


This scenario—the details of which depend somewhat on the boundary conditions (free, 
cylindrical, cyclic, toroidal)—has been expounded in detailed studies of partition function 
zeros [^|9 ,[H ,12 . Indeed, either of the two mechanisms for eigenvalue cancellations outlined 
above have an incidence on the accumulation points of partition function zeros, by the 
Beraha-Kahane-Weiss theorem. In particular, the chromatic line v = —\ intersects the BK 
phase for 0 < Q < 3 on the square lattice, so the accumulation of chromatic zeros is conhned 
to this region. For Q > 3 the chromatic line renormalises to infinite temperature {v = 0) 
and is hence non-critical. 

Remarks. 1. The lower branch of ([^ (i.e., ubk = —VQ) is the analytic continuation of the 
upper branch of (|^, viz., of the critical ferromagnetic curve vp. 

2. The 3-state Potts antiferromagnet on the square lattice is critical on the chromatic 
line V = —1, but when the temperature is non-zero (i.e., —1 < n < 0) it is disordered; so it 
renormalises to n = 0. 


2.3. The Berker-Kadanoff phase 

The BK phase actually exists on any lattice and governs a hnite part of the AF (n < 0 and 
Q > 0) part of the phase diagram for the following reasons, which were succinctly exposed 
in (4^. First, the ferromagnetic transition curve vp —generalising the upper branch of 


must exist on any lattice, by the universality of the CG description with ^ < g < 1. In the 
(Q,'y)-plane this curve must have a vertical tangent at the origin, since the corresponding 
symplectic fermion CFT (with central charge c = —2) is known to describe spanning 
Indeed, if that tangent was not vertical, we would have instead a spanning- 


trees 


43 


forest model 44 , which corresponds to perturbing the free fermion CFT by a four-fermion 


term 45 , which would render it non-critical 45,46 


This means that the critical curve Q{v) has a vanishing derivative Q'(0) = 0 on any 
latticejf] Barring the (unlikely) accident that also Q"{0) = 0, we infer that the critical curve 


f This can indeed be checked explicitly for the square and triangular lattices, where the exact critical curve 
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continues analytically from the first quadrant Q,v > 0 into the fourth quadrant Q > 0, 
v < 0. Invoking again the universal CG description, now with g < we have xr > 2 by 
0. and a hnite range of x-values will indeed be controlled by the BK critical curve ubk- 
This BK phase is bounded in the {Q,v) plane by two smooth curves v± such that for any 
Q € [0,Qc), 0 > v+{Q) > V-{Q). The two curves merge at some value Qc > 0, signaling 
the termination of the BK phase. The inequality Qc < 4 for two-dimensional models follows 


from quantum group results 32 


The upper curve is usually identihed with the critical AF curve, = uaf, and hits the 
point (Q, v) = (0, 0) at a hnite negative slope. This value is related to the critical coupling 


of the corresponding spanning-forest model 44 


Some unusual features of the BK phase have already been pointed out in the previous 
section about the square-lattice model, but they hold true for any two-dimensional lattice: 
At Q = Bk, there are vanishing amplitudes and eigenvalue cancellations, so that the actual 
ground state is deeply buried in the spectrum of the corresponding TM. We now present 
another related argument, this time directly in the continuum limit, leading to the conclusion 
that the BK phase should be dismissed as unphysical. The critical behaviour of the BK 
phase can be obtained by analytic continuation of the CG results for the ferromagnetic 
transition {g > ^). It turns out 47 that some of the critical exponents (namely, magnetic 
and watermelon exponents—to be dehned in the next section) become negative when g < 
This means that the analytic continuation of the ferromagnetic ground state is not anymore 
the lowest-energy state in the BK phase, or alternatively, that correlation functions will 
not decay but rather grow with distance. Such features are clearly unphysical and non- 
probabilistic. 


The existence of the BK phase has been verihed for all Archimedean lattices in 14,15 


and its extent in the (Q,n) plane has been accurately estimated. It was found that in most, 
but not all 13 , cases Qc = 4. 

The critical behaviour on the chromatic line v = —1 (for general Q G M) then depends in 
a crucial way on its intersection with the BK phase. We have seen above that this intersection 
is the interval 0 < Q < 3 for the square lattice. The remainder of the interval [0,4], namely 
3 < Q < 4 in this case, will then have a different behaviour. For the square lattice this 
behaviour is just non-critical. The triangular lattice however offers a much more interesting 
alternative, as we shall see presently. 

Remark. For non-planar recursive lattices [^, the BK phase is also conjectuted to exist, 
and the value for Qc might be greater than 4. 


2.4- The triangular lattice 

The {Q,v) phase diagram on the triangular lattice appears to be considerably more 
complicated than was the case on the square lattice. In particular the part close to 


is known; see Q and (12). 
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Triangular-lattice Potts model 



0 12 3 4 
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Figure 2: Phase diagram of the triangular-lattice Q-state Potts model in the real {Q, v) plane. 
The cubic (12) corresponds to the three thick lines labeled up, '^^bk, and n_, respectively. 
The horizontal line corresponds to n = —1 (13). The curve between up and ubk is our 
best numerical estimate of the AF curve = uaf for this lattice (4^. This curve has a 


bifurcation point T, leading to two branches: one goes to the point {Qo, —1) [cf. (19)], while 
the other goes to the critical point (4, —1). The region dehned by these three points and 
(4, —1.253) is governed by the regime IV that we will consider hereafter, while the BK phase 
corresponds to regime I. (Compare to figure]^) The vertical dashed lines correspond to the 
Beraha numbers [cf. (|^], shown here for 3 < fc < 11. The solid dots depict the known 
critical points for integer values of Q. 


Q = A possesses several intricate features, as witnessed by numerical studies of different 

This model is integrable along the three branches of the cubic [4^ 


types 10-12,44,48 


-|- 3n^ — Q . 


( 12 ) 


Fortunately, the chromatic line 


V = —1. 


(13) 


is also integrable 17,18 and understanding it in detail should provide valuable information 
that will also be useful for disentangling the full {Q,v) phase diagram. 

A schematic phase diagram which is compatible with our current understanding is shown 
in hgurej^ Let ns parameterise the cubic curve (12) by 

'27r(l -^)' 


V = -1-1-2 cos 


Q = A cos"^ {7rg) 


(14) 


where — | < S' < 1. The range 0 < s < 1 then has then same CG interpretation as 
for the square lattice: ^ < g < I corresponds to up, and 0 < s < ^ to ubk- The interval 
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— I < g < 0 (corresponding to n_) cannot be interpreted within this CG, but has been shown 
numerically to belong to the same universality class 38p^ as the antiferromagnetic curve 
(|^ for the square-lattice Potts model (with the same value of Q). 

This implies that the middle branch ubk of (12) will control the BK phase, which 
should extend from the lower branch of (12), that we hence identify with n_, up to some 
AF transition curve uaf = that unfortunately has not yet been determined analytically. 
This curve cannot be the chromatic line (fl^, since the latter is not always above the middle 


branch of (12). At the origin, {Q,v) = (0,0), the AF curve describes a spanning-forest 


problem and its slope is known numerically . It has be followed numerically to higher Q 
in [^, and other pieces of yet unpublished work, and is depicted in £gure|^ 

Along the chromatic line v = —1, Baxter has found three distinct analytic 
expressions for the TM eigenvalue gi, whose logarithm is the bulk free energy (per site): 


9i{Q) 

log^2(Q) 


log^3(Q) 


1 (1 - x6^-3)(l - a;6^-2)2(i - 

X (1 — x®'5“®)(l — x®'?“"^)(l — a;®'?)(l — 


,, sinh kO 
ak ;— 

-no 2A; 


dk 


sinh[fc(7r — 29)/2] 

(2 cosh kO — 1) sinh(7r/c/2) 

cosh[/c(7r — 29)/2] 

(2 cosh k9 + 1) cosh(7rfc/2) 

sinh k9 [sinh k{n — 9)] 


k sinh Tik [2 cosh A;(7r — 6*) — 1] ’ 


where we have parameterised 


Q = 2 — X — X ^ = 2 + 2 cos 9 , 


(15) 


(16) 

(17) 


(18) 


with |x| < 1 and 0 < Re(6*) < tt. (See Ref. 17 for details.) 

In his second paper [^, Baxter determined the regions of Q € C for which each of the 
expressions (15)-(p(^ is dominant (see figure [^. This was done by investigating the ratios 
g 2 /gi and g^/gi, which can both the expressed as inhnite products, by a combination of 
symmetry arguments and numerical analysis. His hnal result is: 


(i) gi is dominant for Q G (—oo, 0) U (4, oo). 

(ii) g 2 is dominant for Q G (Qo,4). 

(iii) 5 f 3 is dominant for Q G (0,Qo)- 


where 


10,18 


Qo = 3.819 671731239 719-■■ . 


(19) 


Returning now to the phase diagram in the (Q, n)-plane, the simplest scenario that can 
account for these findings is that the AF transition curve must split before Qq into a branch 
that goes down to (Q,n) = (Qo 5 ~l) and delimits the BK phase, and another branch that 
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extends further to the 4-colouring point {Q,v) = (4,-1). This scenario, which is shown 
schematically in hgure is conhrmed by studies of the limiting curves and the critical 
polynomial [^. The region n < — 1 appears to contain further intricate features, which are 
beyond the scope of the present discussion. 

The region (i) where gi dominates is non-critical. In that region, exact expressions for 
the surface and corner free energies are also known 29 


of angle | and ^ can even be determined separately 


30 


and the contributions from corners 
. In the limit Q —)■ 4+, these corner 


free energies develop essential singularities that have been determined from the asymptotic 
analysis of the exact expressions |^ . 

For the remainder of this paper we shall restrict the attention to the critical regions (ii) 
and (iii). For reasons that will become clear from the mappings in section]^ to a spin-one 
vertex model, we shall refer to the critical theory that dominates in region (iii) as regime I, 
and to the theory dominating (ii) as regime IV. (See hgure [^) Our objective is to extend 
Baxter’s analysis to compute the exact critical exponents in either of these regimes, and 
determine the corresponding CFT. 

Regime I is relatively straightforward to deal with, using the mapping to the 0(n) 
model 19 and the BA analysis of the associated 7-vertex model 21 . Choosing a particular 


twist of the vertex model to ensure the equivalence with the Potts model will enable us to 
compute the corresponding critical exponents; see section 4^ The end result is that regime I 
is also described by the BK universality class. It is remarkable that the triangular-lattice 
Potts model has two distinct integrable curves, (13) and the middle branch of (12), with the 
same critical behaviour. 

Regime IV is our main interest here. It cannot be obtained within the 7-vertex 
model 21 , and indeed in that reference the existence of the region Q G (Qo,4) was relegated 
to a terse footnote. However, using further mappings to a 19-vertex model (see section]^, we 
shall obtain an efficient handle on regime IV and derive the corresponding critical exponents 
and CFT interpretation. Achieving this is going to be the principal motivation for the 
remainder of this paper. 

Remarks. 1. Although Baxter’s analysis of the eigenvalues ([T^-([T^ has been subject of 
some debate —in particular concerning the precise numerical evaluation of the inhnite 
products—his conclusions for Q G M have withstood a number of stringent numerical 


tests . 

2. In the square-lattice case, we could say that Qq = 3. The corresponding chromatic 
polynomial does not dehne an integrable theory. However, since the region 0 < Q < Qo is 
attracted to ubk under the renormalisation group flow, it can be identihed with regime I. 
The behaviour right aX. Q = Qq is given instead by the AF curve n+. 

3. The extension of Baxter’s work to the computation of critical exponents, by means of 
a BA analysis, is deferred to section]^ where a twist should be included in the BA equations. 
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Baxter 17 found an ingenious coordinate BA for the triangular-lattice chromatic polynomial, 


involving several kinds of particle trajectories (that he represented as dotted, full, double-full 
and wiggly lines). In view of the subsequent developments in the classihcation of integrable 
systems, it appears preferable to instead bring the chromatic polynomial into contact with 
a (by now) well-studied integrable model to which we can apply all the (by now) standard 
BA machinery. In this section we therefore reexpress the colouring problem in terms of an 
integrable, spin-one vertex model on the square lattice, known as the Izergin-Korepin (IK) 
(or a^ 2 ^) model. 

The route that we shall take passes through quite a number of different formulations 
as intermediate steps. We emphasise that although each of the mappings exhibited here 
is exact, the various models appearing have different domains of validity. Some require 
Q G N, while other allow for any Q G M; sometimes the twist (that determines the weight of 
non-contractible clusters and loops) is a free parameter, and sometimes it is hxed. In more 
technical terms, each model provides a distinct representation of the underlying Temperley- 
Lieb (or dilute TL) algebra, which is not always faithful with respect to the most general 
formulation (allowing for arbitrary values of Q and the twist). It follows that some of the 
eigenstates—and hence some of the critical exponents—are absent from some of the models 
presented, and not all of the correlation functions can be directly interpreted in terms of the 
original colouring problem. To keep track of these subtleties, our discussion puts a special 
emphasis on the dehnition of the corresponding transfer matrices, as well as on the various 
physical observables of interest. 

Another key point is to keep carefully track of the twist. The bulk free energy obtained 
by Baxter is independent on the twist, so he did not need to introduce it. However, since 
we are eventually aiming at exact expressions for the critical exponents—which are sensitive 
to the boundary condition dehned by the twist—we need to ensure that each model in the 
series of mappings is identically twisted. 

It is worth stressing that we will start from the colouring problem dehned on a strip 
of the triangular lattice with cylindrical boundary conditions and a width L multiple of 
3. This choice allows for a 3-colouring for such strips, with vertices of each of the three 
sublattices having constant colour. Strip graphs of widths L not a multiple of 3 correspond 
to a topologically frustrated ground state, depending on the value of L (mod 3), and will be 
considered elsewhere (SOl. 


3.1. Colouring problem 

We start from the Potts model dehned on a strip of the triangular lattice G = (V, E), whose 
partition function we recall to be given by 




K&a a 
e *’ j 


Wi} {i,j} 


( 20 ) 
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where Q € N and G M. The colouring problem is recovered for K ^ —oo. Viewing the 
longitudinal direction of the strip as imaginary time, we can build the system row by row 
acting with the TM, which we represent as follows: 



where the dashed lines indicate that the two ends are to be joined to recover the cylinder 
geometry. 

For generic values of K, the TM in this spin representation acts on the space of all 
spin conhgurations of a given row, {1,2,... ,Q}^. Therefore, its dimension is given simply 
by Q^. In the colouring problem (K = —oo) this dimension is reduced because of the 
nearest-neighbour constraint dj 7 ^ aj. In terms of the corresponding colouring matrix 
A = (1—the dimension reads di = Tr A^. Now A has one eigenvector Vi{ai) = 1 
with eigenvalue Q — I, and Q — I eigenvectors Vkici) = for A; = 2, 3,..., Q, each 

with eigenvalue —1. Therefore di = {Q — 1)^ + {Q — which obviously grows 

asymptotically as {Q — 1)^. Note that di can be further reduced by taking into account the 
(transverse) translational symmetry of the strip. 

The main drawback of the TM method in the spin representation is that for each value 
of Q G N we have to produce a distinct TM. 

We now turn to the physical correlation functions and critical exponents of interest 
in the colouring problem. The first of these is the so-called magnetic correlation function, 
dehned between two sites i,j of the triangular lattice as 



On the critical lines of the phase diagram, the magnetic correlation function has an algebraic 
asymptotic behaviour: 


~ \i-j\ for|i-jj>l, ( 22 ) 

where xh is the magnetic exponent. Away from critical points, it decays exponentially fast 
with the distance, GH{i,j) ~ with a typical length scale given by the correlation 

length The asymptotic behaviour of the correlation length in the vicinity of a critical 
point allows for the dehnition of another critical exponent: 


rsj 


V — V, 


— V 


(23) 
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which is related to the thermal exponent xt through 

1 


10 = 


2-.r- '''> 

Two-point functions of more general operators, acting on groups of more than one spin, 

They include in particular the so-called /c-cluster ‘watermelon’ 


have been defined in 51 


operators, which can however be formulated more easily in the FK representation to which 
we turn next. 


3.2. FK representation 

The partition function (20) can be rewritten in the FK representation as follows 22 : 

ZS^lQ.v) = ^ , (25) 

ACE 

where n = e^ — 1 = —lin the colouring problem, and k{A) is the number of connected 
components of the spanning subgraph {V,A). As this function is a polynomial in both Q 
and V, we can extend those variables analytically to arbitrary real or even complex numbers. 
The original colouring model is in principle recovered when we specialise (25) to n = — 1 and 


Q G N; however one must ensure that the free energy and the correlation functions obtained 
in the cluster representation still make sense in the colouring problem. Whether this is true 
or not depends, as we will see, on the regime under consideration. 

As in the spin representation, we can define the FK representation of the TM on the 
triangular lattice, with the important difference that the latter considers all values of Q G C 
at the same time. The computation of the TM for a triangular-lattice strip with cylindrical 
boundary conditions is not straightforward: one needs to consider first a strip of width L -|-1 
with free transverse boundary conditions, and at the end, identify sites 1 and L-l-1 [8 18 . The 
final TM acts on the space of all non-crossing, non-nearest-neighbour partitions of the set 
{1, 2,..., L} on a circle. The non-crossing property is due to the planarity of the graph, and 
the non-nearest-neighbor property is due to n = —1 (i.e., neighbour spins cannot be coloured 
alike) [8,29 . The number of such partitions is given by the Riordan number for L > 2; 
its asymptotic behaviour is Rl = 3^L“^/^3a/3/(8i/F)[ 1-|-0(1/L)] (see e.g. ref. [^ section 3.3 
and references therein]). If we also take into account the translational symmetry along the 
transverse direction, the number of partition classes modulo translations is asymptotically 
equal to Rl/L. 

Turning to physical observables, the magnetic correlation function (21) gets 
straightforwardly translated in the FK representation as 


GA'iC) = 


G 


Q Z, 

0-1 zp-{Q,v) 


(26) 


where is the partition sum of FK configurations where sites i and f belong to the 

same cluster. In the TM picture, the points i and j get mapped onto the infinite top and 
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infinite bottom of the cylinder, and the exponent xh can be obtained from the scaling of the 
TM eigenvalues as 


log 


A, 


{L) 


A 


{L) 


2tiGxh (1 


(27) 


where G = ^ is a geometrical factor 10 , Ag^^ is the largest TM eigenvalue, and A^^^ is 


the largest eigenvalue in the subsector where one cluster (considered ‘marked’) is imposed 
to propagate along the vertical direction of the cylinder, and hence cannot be destroyed by 
the action of the TM. This dehnition can actually be extended, as we can enlarge our state 
space by imposing the propagation of k distinct marked clusters from the bottom to the top 
of the cylinder, leading to the set of exponents 


log 


A, 


(i) 




271:0 Xk (1 

—u 


(28) 


The case k = 1 actually contains two distinct situations, depending whether the propagating 
cluster is allowed to wrap around the cylinder or not; we will come back to this issue in 
section 13.41 

In a similar fashion, we expect that the thermal exponent should be related to the 
relative scaling of the subdominant eigenvalue Ag^^ in the /c = 0 sector, with respect to the 
dominant one Aq^\ namely 


log 


A, 


{L) 


a; 


/(L) 


27:0 xt (1 


(29) 


Remark. It is worth mentioning that for Q = Bp [cf. (|^], with integer p > 3, the 
Potts-model partition function allows for another probabilistic interpretation (with positive 
weights), namely in terms of a restricted-solid-on-solid (RSOS) model 37 52 . We will discuss 
this issue in section 1 


3.3. Duality transformation 


Given a Q-state Potts model with parameter v in the FK representation (25) and dehned 


on a triangular lattice O = {V,E), a duality transformation allows us to obtain an 
equivalent Q-state Potts model on the dual, hexagonal lattice O* = (y*,E*), by putting a 
FK link on each edge of E* not crossed by a FK link on E (see hgurej^. 

The corresponding dual parameter is given by 


= - = -Q. 

V 


( 30 ) 


and we can represent the corresponding TM as 



















Q-colourings of the triangular lattice 


17 



Figure 3: Configuration of FK clusters on the triangular lattice G (black), and the associated 
conhguration of dual FK clusters on the hexagonal lattice G* (blue). In red, we have 
represented the surrounding loops on the medial, kagome lattice. 



where the dual, hexagonal lattice has been represented in blue. The space of states in terms 
of a set of connectivities follows from a construction analogous to that used on the triangular 
lattice, and we point out that imposing a number k of propagating clusters in the triangular- 
lattice Potts model is equivalent to imposing the same number of propagating clusters in 
the dual, hexagonal-lattice model, so that the dehnition of critical exponents sketched in the 
previous section translates straightforwardly to the dual model. 

Note that under duality the chromatic polynomial xg of a planar graph G transforms 
into the flow polynomial on the dual graph G*. To be more precise, Xg*{.Q) = Q ^g{Q)- 

3.4- Six-vertex model on the kagome lattice 

The Potts model on the dual, hexagonal lattice G*, can be further mapped mapped onto 
a model of loops living on the edges on the medial, kagome lattice, and corresponding to 
the surrounding contours of the FK clusters in the spanning subgraph (K*, A) (see £gure[^. 
This is achieved through the following identity 

Zl^iQy) = 

AGE* 


(31) 
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X 


* 



(32) 


The conhgurations of this loop model can now be mapped onto those of a six-vertex 
model on the kagome lattice 49 . This is done by reexpressing the factor terms of 


local factors by setting [cf. (18) 


Q 

•/Q = 2 cos-, 


(33) 


and summing over the orientations of loops, as represented in the hrst two lines of hgure 
Notice that in this hgure we represent only one of the possible three types of vertices on 
the kagome lattice (related by 27r/3 rotations); however, contrary to what happens for the 
square lattice, we can always see which weights are which because the two angles in the 
kagome lattice are distinct: tt/S and 27r/3, respectively. 

The TM of the resulting six-vertex model can be represented as follows (the vertices 
and edges of the medial kagome lattice are depicted as red dots and lines, respectively): 



The degrees of freedom are now the orientations on each edge of the kagome lattice, so 
the TM acts on each row on the space of 2L edges oriented either upwards and downwards, 
therefore with a total dimension 2^^ = 4^. 

As suggested by the hgure above, a twist factor e^*^ has to be introduced in the six- 
vertex model in order to ensure the correspondence with the original periodic Potts model. 
It is represented by a red wavy line, which contributes a factor e ^®2 each time it is crossed 
by a rightwards (resp. leftwards) arrow of the six-vertex model. The twist parameter ip must 
be chosen as follows: 

• In the sector with no propagating clusters (/c = 0), there may be closed loops on the 
kagome lattice, that wind around the transverse periodic boundary condition. To recover 
the original model, these loops must be given the same weight ^/Q = 2 cos | as the one 
given to contractible loops. Choosing the twist f ~ f shall see, | 

depending on the regime), since each loop has to be summed over the two possible 
orientations, it is given a weight 
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as required. Therefore we expect the central charge of the Potts model to be given in 


terms of the free energy (density per site) by 55,56 


fM = ^ log 




ttGc (I 

- /bulk + -^+ol ^ 


(35) 


where A.^Q\ip) is the largest eigenvalue of the six-vertex TM, appropriately twisted, in 
the sector with £ = 0 ‘through-lines’ (Potts loops propagating along the vertical direction 
of the cylinder). 

In the sector with a single propagating FK cluster {k = 1) we have two distinct 
situations: 

— If the propagating FK cluster is allowed to wrap around the cylinder, then this 
translates in the six-vertex language to no loop propagating in the vertical direction, 
but no closed loop allowed to wrap around the cylinder. Forbidding loops on the 
kagome lattice to wrap around the cylinder amounts to giving such loops a weight 
0, therefore to fixing the twist to ^ = |. The magnetic exponent can thus be 
obtained as 


log 




K 




2710 xh f 1 
L L 


(36) 


— If the propagating FK cluster is not allowed to wrap around the cylinder, then 
it translates in the Potts loops/six-vertex language into the presence of two loops 
propagating in the vertical direction (‘through-lines’) delimiting the vertical borders 
of the propagating FK cluster. The situation is then similar to the case k > 2, 
described in the following paragraph. 

In the sectors with k > 2 propagating clusters (or k = 1 cluster not allowed to wrap 
around the cylinder, see above), there are no non-contractible loops anymore. Each 
propagating cluster is associated with two Potts through-lines, and one can take each of 
the i = 2k through-lines around the axis of the cylinder without affecting the Boltzmann 
weights, if every such line acquires in doing so a phase which is a root of unity. This 
means that the eigenvalues of the original model are obtained in the twisted vertex 
model with a twist being of the form with any integer p satisfying gcd(p, £) = 1. 


3.5. 0{n) model on the hexagonal lattice 

As shown in the last line of hgure the configurations of the six-vertex model can in turn 
be mapped onto those of dilute loops living on the edges on the hexagonal lattice. In this 
mapping, the thin spin-| loops run around the inside and outside of the thick spin-1 loops, 
therefore a dilute loop on the hexagonal lattice is a double Potts loop on the kagome lattice. 
The resulting gas of dilute loops is nothing but the 0(n) model on the hexagonal 


lattice 19 , where the weight of each closed loop is 

n = —2 cos 20 


(37) 
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Figure 4: 

Mapping from the Potts model (31) on 

the hexagonal lattice 

(hrst line) onto a 


six-vertex model on the kagome lattice (seconc 
on the hexagonal lattice (third line). 


line), and then onto a model of dilute loops 


and the fugacity of each loop segment is 

K 


1 

2sinf 


The corresponding TM can be represented as follows: 


(38) 



One needs to be careful about what becomes in this model of the twist parameter 
introduced in the six-vertex model: each loop segment of the 0(n) model crossing the 
vertical seam has to be summed over two possible orientations. Now, from hgure one 
sees that each oriented 0(n) loop crossing the vertical seam is equivalent to two edges of the 
six-vertex model crossing the seam with the same orientation, hence a twist ip rather than 
^ (where p has to be conveniently chosen depending on the sector we are considering, as 
explained in the previous section). 

It is important to stress that making the right choice for the twist in order to ensure 
the correspondence with the original periodic Potts model results in a ‘wrong’ choice of the 
weight for non-contractible loops in the 0(n) model. For instance, in the k = 0 sector, if we 
choose I = |, we get the right weight for the winding FK clusters ^/Q = 2cos|; but the 
wrong one for the winding 0(n) loops n = 2cos6'. Similarly, if we choose ^ we get 

the right weight for a winding 0{n) loop n = —2 cos20, but the wrong one for a winding 
FK cluster = 2sin0. In other words, the periodic Potts model is equivalent to a 0(n) 
model on the hexagonal lattice with a particular choice of twisted boundary conditions. 
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Pi P2 P3 Pi Pb Pe P7 Ps P9 



Figure 5: Vertices of the square lattice 0(n) loop model (first row), together with the 
corresponding Boltzmann weights, and their mapping onto the vertices of the hexagonal 
lattice 0(n) model when the weight pg is zero. 


The TM acts on the set of non-crossing dilute matchings of the set {1, 2,..., L}. By 
dilute (in opposition to perfect) we here mean that any given point may be empty and hence 
not paired with another point. The number of such matchings is given by the Motzkin 
number Ml for L > 1; its asymptotic behaviour is 3^L“^/^3 a/3/(2a/F)[1 -|- 0{1/L)\. 


3.6. Spin-one vertex model on the sguare lattice 


The 0(n) model on the hexagonal lattice is a particular case of the 0(n) model on the 
square lattice, whose vertex weights are parameterised as shown in the hrst line of hgure 
Indeed, in the particular case pg = 0, the vertices can be stretched horizontally by adding 
one additional horizontal edge as shown in the second line of the hgure. 

An integrable choice of the Boltzmann weights is well known (we use here the 
notations of (^): 


Pi(w) 

P2{u) 

Pi{u) 

Peiu) 

Psiu) 

Pdiu) 


1 + 


2 cos(n -|- y) sin u sin | 
(2 cosP — 1) sin^ 9 


P3{u) = 
P5{u) = 
Priu) = 


COs(u -h y) 
COS y 

sinu 


cosf ’ 


2 cos(m -|- y) sin M sin | 
(2 cos 9 — 1) sin^ 9 


sin(2M + f) 


sm 


0 


2(cos 9 -|- cos 29) 
0 


cos m 


sinu 


2 cos^ I sin I (2 cos 9 — 1) 


(39a) 

(39b) 

(39c) 

(39d) 

(39e) 

(39f) 


In all the following we rescale p^ and ps by a factor —1, which is just a change of basis. This 
model has several symmetries: 
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• The transformation u ^ ^ ^ — u for fixed 6 is equivalent to rotating the vertices of 

hgure through | (crossing symmetry). 

• Up to some gauge symmetry, the weights are also invariant under vr shifts of u, as well 
as under 27r shifts of 9. 


For a given 9, there is only one non-trivial value of the spectral parameter u such that 
namely u = f — f • The corresponding weights, which we rescale such that pi = 1, 
following [cf. (39) 


P9 = 0 , 
are the 


Pi 

= 1, 

(40a) 

P2 

1 

~ ~ r, ■ 6» ’ 

2sm 1 

(40b) 

P4 

1 

— ■■■ — Ps — A ■ 2 e 1 

4 sm 1 

(40c) 

P9 

= 0, 

(40d) 


so they recover precisely the fugacity K given by eq. (38). 

Now, it is well-known (see 58 for a review) that the integrable 0(n) model on the square 
lattice can be reformulated as a spin-one vertex model, the so-called IK (or a^^^) model. The 
equivalence is worked out by once again giving orientations to the loops (to be summed 
over), and turning the loop weight n into local, angular contributions. Each occupied edge is 
associated with a state ±1 (depending on the orientation of the corresponding loop), while 
empty edges are associated with the state 0. 

We can now bend the lattice, and represent the corresponding TM as follows: 



where each edge carries a = —1, 0,1 state, the interactions are encoded by local vertex 
weights depending on the 19 possible states of the adjacent edges, and the vertical seam (red 
wavy line) corresponds to a boundary twist 1,0**^. The fact that the vertical seam now 
only crosses one horizontal edge of each row (whereas it also crosses the last “vertical” edge 
in the six-vertex formulation on the kagome lattice) is irrelevant, since the vertical crossing 
can be eliminated by a change of basis. 

The TM of the model commutes with the total magnetization 

i=l 


m 


(41) 
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which can naturally be associated to the number of through-lines in the 0(n) loop 
formulationjl] Therefore each sector of fixed m can be considered independently. 

Rather than checking now explicitly that the eigenvalues dominating the spectrum of 
the TM for the different sectors of the colouring problem are indeed recovered by the spin- 
one vertex model, we defer this discussion to the next section, where we will deal with the 
BA construction of these different eigenvalues. 

The total dimension of the spin-one vertex model TM, summed over all magnetisation 
sectors m, is obviously 3^. In the m = 0 sector the dimension is given by the central 
trinomial coefficient, that is, the largest coefficient in the expansion of (1 -|- a: -|- This 

has asymptotic behaviour 3'^L“^/^'^3/47r[l -|- 0(1/L)]. Similarly, the dimension in sector 
m = 1 is given by the next-to-central trinomial coefficient, and more generally, that of the 
sector with magnetisation m by the trinomial coefficient that is m places away from the 
centre. 

We stress that these dimensions are larger than those of the 0(n) loop model, for 
the reason that the TM of the latter does not distinguish between contractible and non- 
contractible loops (it gives a weight n to either type of loop). It is however possible to work 
with an ‘augmented’ loop model TM in which each arc connecting two points is distinguished 
by the parity (even or odd) of the number of times it has crossed the vertical seam. When 
closing an even (resp. odd) arc we obtain a contractible (resp. non-contractible) loop. It 
is easy to see that states of this augmented loop model are in bijection with those of the 
spin-one vertex model. It follows by the above equivalences that the spectra of the respective 
TM are identical. 


4. Bethe Ansatz solution 


After all the mappings of section we have managed to reformulate the original colouring 
problem as an integrable spin-one model, namely the ^ model with a particular value of 
the twist, related to the weight of non-contractible FK clusters. This means that we are now 
in a good position to apply the toolbox of the BA method. 

The goal of the present section is to set up the general framework, and define properly 
the regimes I and IV that were introduced in section 2.4[ as well as the possible excitations of 
the ground state configuration at finite size L. With these ingredients at hand, we can then 
in section take the continuum limit and work out the corresponding critical exponents. 


f.l. Bethe Ansatz eguations 

It is well-known how to construct the eigenstates of the integrable spin-one model through 
the BA. Each eigenstate is parameterised by a set of Bethe roots, solution of the Bethe 
Ansatz equations (BAE). In our particular case, the BAE for the square-lattice 0(n) model 

f Here and below, the word through-line refers to the spin-one 0(n) loop model. Each such through-line 
corresponds to a pair of through-lines in the spin-1/2 loop model; see figure bl 
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(which is related to the IK model) are given by: 


f sinh {Xj - if) 
Vsinh(Aj +if) 


= e^^ 


n 


sinh {Xj — Aj — i^) cosh (Aj — Aj + i|) 




'r sinh {Xj — Aj + i^) cosh (Aj — Aj — i|) 


(42) 


From there, we look for the leading TM eigenstates of the colouring problem in terms 
of Bethe roots for an appropriate twist. The following conclusions have been known since 
the paper of Baxter 


The range 0 < Q < Qo is described by the so-called regime I in the notations of 59 
Namely the ground state is made of roots with imaginary part 
The Qo < Q < 4 regime is dominated by another conhguration of roots 


17 


, namely ^ 

real roots, and ^ roots with imaginary part |. It is worth noting that the corresponding 
eigenstates are absent from the loop formulations, both on the square and hexagonal 
lattice. Also, it never dominates at one of the isotropic values of the spectral parameter 


studied in 58,59 , nor does it dominate the spectrum of the corresponding Hamiltonians, 


obtained by differentiating the TM with respect to the spectral parameter u, at the value 
u = 0. As a consequence, this regime has escaped the attention of previous works on 
the IK vertex model, in which three regimes labelled I, II and III were identihed and 
studied 58,59 . Extending this nomenclature, we shall call the regime dominated by 


this root conhguration regime IV. 


Before going further in the study of the regimes, let us compare the BAE (42) with 


those of 17,21 , where a 7-vertex formulation of the hexagonal-lattice 0(n) loop model is 
solved without any reference to the IK 19-vertex model, which, as we have seen, is related 
to the more general square-lattice 0(n) loop model. The BAE are written in terms of a set 
of complex Bethe roots {zj}, and read 

S{zi,Zj) 


= (- 1 ) 




n 

h¥^j) 


S(zj,Zi 


(43) 


where the scattering phases are 

S{z,w) = (l — z — w + zw + t^z) {l — 2w + zw + t^w) , 
and t^ = 2 — 2cos9. In (43), we use the parameter n = L + — n+ equal to the number 


(44) 


of lines in Baxter’s formulation 17 , with n+ (resp. n_) being the number of upward (resp. 
downward) pointing vertical arrows in each row 1^. Our sector k = 0 corresponds to 


these authors 


21 


?T,_ = n+, hence n = L. Following |21|, we set Zj = where r = e*®; however, while 


further dehne Wj 


^ r+Wj 

, we choose instead to set Wj = 


,2A, 


From there, 

the BAE (43) are easily seen to take the form (42), where the correspondence between the 


twist parameters is simply cp = e. In other words, the BAE underlying the construction of 
eigenstates for the two models (the 7-vertex model of and the IK 19-vertex model) are 
the same, even though the space of states of the two models are not equivalent j|] 

f In particular, we expect that either model is able to reproduce both regimes I and IV. Indeed, the BAE 
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Regime I, originally investigated in the context of the square lattice [^, was studied in great 
detail for the 0(?7,) loop model on the hexagonal lattice in 21 . The CFT is that of a single 
compactihed bosonic held, and in particular, the dependence of the central charge c on the 
twist (f was exactly found to b<i 


c((f) = 1- 


n6 


(45) 


while the set of exponents with respect to the untwisted central charge c = 1 were found as 

m‘^9 


Ar 


,vr 


to +"’»• 


and therefore, with respect to the twisted theory (45), 

m‘^9 


.TT 


^m,w — 




47r 


+ 9 "toO’ 


(46) 


(47) 


where m the total magnetization introduced in (41), and w can be thought of as related to 
the electric charge in the CG picture. 

It turns out that the ground state of the 0(n) model with twist ip = 9 is not observed 
in the Potts-model spectrum, and the central charge of the Potts model is recovered by 
choosing the twist in this regime as (p = 27r — 6*jjj] Then ( [4^ reduces to 


CPotts 1 


3{27i-9f 


= 1 - 


6(p — 1)^ 


(48) 


7i9 p 

where we have parameterised 9 = 2tt/p. This is precisely the central charge for the BK 
phase found by Saleur 
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47 


, matching the CG results for coupling ^ G (0, |). The 
Gleg watermelon exponents are straightforwardly obtained from (|47j), namely 


X e Q 
2 


f9 {271 - 9y 


{p-lf 


(49) 


IbTT 4:Ti 9 8p 2p 

and coincide with the known values for bulk watermelon exponents in the BK phase (see for 
instance 60,^). As explained in section 3.4, the magnetic exponent is associated to the 


scaling of the TM eigenvalue with twist | = |, namely 


(*^Potts 


C(^)) = ^ 


p {p- ly 


2p 


(50) 


construction from the 7-vertex model is the one originally used by Baxter [17| , which gave the first evidence 
for the existence of regime IV. 

§ Note that there is a typo in eq. (2.62) of 


21 


. We give the right expression here. 


II We have already remarked (in section 3.41 that both possibilities for the twist parameter would give the 
correct weight to non-contractible loops, but the comparison with the spectrum of the Potts model here 
imposes the latter choice as the correct one. 
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which once again coincides with the known value in the BK phase. The thermal exponent is 
obtained from the scaling of the subdominant eigenvalues in the ground state {m = 0,w = 0) 
sector. These are obtained by taking values of the twist parameter shifted by 27r, namely 
(f = 27r(l + j) — 9, resulting in 

= A (cpotts - c(27r(l + j) -9)) = J [{j + 2)p - 2] . (51) 


The two most relevant of these exponents (excluding the case j = 
spectrum) correspond to j = 1 and j = —2, hence respectively 


( 1 ) 

- Jbrp 


= ^-1 > 2 , 

2 - ’ 




(- 2 ) 


Recalling that P = the first of these coincides with eq. (11). 


1, not observed in the 

2. (52) 


4 . 3 . Regime IV: structure of the leading eigenlevels 


We now study in detail the so-called regime IV. Let us set notations: For sizes L G 3N, the 
ground state is described by mi = ^ real roots Xi, and m 2 = f roots with imaginary part 
I that we denote 2 /a: + if with G M. 

The Bethe equations (42) can be recast as: 

L 


( sinh {xj - if) 

Vsinh + if) 


cosh (i/fc - if) 
cosh (l/fc + if) 


= e‘^ 


n 


sinh {xj — Xi — i9) cosh (^Xj — Xi + i|) 


= e‘^ 


sinh {xj — Xi + i9) cosh (xj — Xi — i|) 

cosh {xj — Hk — ’i9) sinh {xj — -|- i|) 

^ cosh {xj -yk + i9) sinh [xj - yu - if) ’ 

sinh {yk - yi - i9) cosh [yk - yi + if) 
sinh {yk - yi + i9) cosh [yk - yi - if) 

cosh {yk — Xj — i9) sinh — Xj -|- if) 

cosh {yk — Xj + W) sinh — Xj — if) 


n 

n 


X 


n 


(53) 


(54) 


Introducing the real function 0(A, a) = ^ log , and taking the logarithm of these 

equations, they can be written as a set of coupled non-linear real equations in the variables 
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Vk}- 


Lcj) y ^j 


( ^ 

Xj, 


V ^’2 


2 


- 




0 {Xj - Xi,d) +(p( Xj - Xi, 


n — 9 


+ ^ 0 (xj - l/fc, 0 + I) + 0 f - Vk,-- 


(55) 


TT + 6* 


^0 [vk, 2 


2 


TrJfc + ^ 


0(|/fc - |/Z,6») + 0 h/fc - i/z, 


Z(^fc) >- 


n — 9 


+ Y1 (l^{yk-Xj,9+^'^ +(j)\^^yk-Xj,-- 


(56) 


where the Jj and Jk are distinct integers (or half-integers), the so-called Bethe integers. In 
the ground state configuration, these are maximally and symmetrically packed, namely, for 
L odd 


-_ _o- - -o- - -0 - - o - - a - -o- Im = ^ 

h _ 1 
6 2 

3 1 I 1 3 LI 


+ 0 1 


_hi 1 ... 

3^2 22 22 32 


Im = 0 


while for L even 


-o- 

--0- 

--0- 

-O- 

+ 1... 

6^2 

1 

1 

L 1 

2 

2 

6 2 

L . 1 3 

3_~'~ 2 _ _2 

1 

1 

3 LI 

2 

2 

2 'G 3- 2 




Im=f 


Im = 0 


From this ground state solution several elementary excitations can be considered: 

• Magnetic excitations, corresponding to removing Smi,Sm 2 roots of each Fermi sea, 
together with a rearranging of the Bethe integers such that these remain symmetric and 
maximally packed. 

• Global backscatterings, corresponding to shifting all the Bethe integers of either Fermi 
sea by an amount 2 wi, or respectively 2 w 2 (where Wi are either integers of half-integers). 

• Particle-hole excitations of either of the Fermi seas, which correspond to shifting 
(resp. (50) of the Bethe roots to a higher/lower Bethe integer. This can be represented 
as follows (for the plus sign): 


o— 
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Other excitations can also be considered, namely roots can be rearranged in complex 
conjngate pairs called 2-strings\ Xi ± icfi. However, snch confignrations can be considered as 
‘artefacts’, in the following sense: not only these conhgnrations are highly nnstable nnder 
changes of 0 or (p and their analytical treatment (if any) escapes the range of the methods 
nsed in this paper, bnt the associated conformal weights can actnally be obtained from the 
nnderstanding of states with no 2-strings. How this comes abont shall be made clearer in 


section |5.1[ where we will examine in detail the strnctnre of the eigenstates of the colonring 
problem in the fc = 0 sector. 


f.f. Regime IV: analysis of the Bethe-Ansatz eguations 


In the continnnm limit, the logarithmic BAE (56) tnrn into eqnations for the densities of 


real parts p{x)^ In Fonrier space (where the densities are denoted by p(a;), ct(c(;)), these 
eqnations are easily rewritten as 

sinh(n (^ - I) _ sinh u (^ - 9) + sinh ^ _ 

^ ' =P+-- -P 


sinh ^ 


sinhf 
sinh ^ 


sinh ^ 


sinhw (f — I) + sinhw^ 


sinh ^ 


a 


(57) 


= a + 


sinhcj (f — 6^) -|- sinh ^ 


sinh ^ 


a 


sinhca (f — |) + sinhca^ 


sinh ^ 


P 


(58) 


These eqnations are easily solved, and we obtain the densities in the thermodynamic 
limit of the real roots in Fonrier space: 


P = 


2 cosh f 
1-1-2 cosh ujO ’ 


a = 


1 + 2 cosha;6* 


(59) 


So, back to real space. 


p{x) = 


1 sinh ^ 


+ sinh ^ 


sinh ^ 


cr 


iy) = 


1 sinh II 
61 a/3 sinh ^ 


(60) 


As discnssed in the previons section, the excitations above the gronnd state are obtained 
by creating ‘holes’ and ‘backscatterings’ in the Fermi seas of x and y (we shall not consider 
here the excitation obtained by forming 2-strings). Introdncing the densities of holes 
and (which play an analogons role as the densities of roots p and a), we can write the 
scattering eqnations in a matrix form 


p + p^ 

a + a‘ 


h 


= s + K 


P 

a 


(61) 
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p=7, 6m^=6'm2=w^=W2=0 



Figure 6: Estimates of the central charge for the ground state (with root conhguration given 
by (72)) for p = 7 (Q ~ 3.2469796) as a function of the twist ip. We show the numerical 
solutions of the BAE for sizes L = 6 (black □), L = 12 (red o), L = 18 (navy blue A), 
L = 24 (pink 0), F = 30 (orange x), L = 36 (violet ♦), L = 42 (dark gray B), and L = 48 
(blue •). The convergence towards the analytic formula c = 2 —^ [cf. (66)], depicted 

as a black curve, is excellent. 


The conformal weights of excitations indexed by the holes numbers { 6 mi, 6 m 2 ) and 


backscatterings {wi,W 2 ) are obtained (at zero twist) as 62 


Ao + A„ = 6m2) {I - K) 


6mi 

61712 


+ {wi, W2) (1 - K) 



(62) 


where the label a denotes for the set of quantum numbers of a given excitation, and in 1 — 77 
the zero-frequency limit a; —)■ 0 has to be taken, so 


l-K 




(63) 


Therefore, 


Aq -|- Aq, 


1 

2 



{5mi^ -I- 5mi^) 



6 mi6m2 


- (wi^ + W1W2 + H- 7 

3 ^ ^ p — 4 


(Wl -7 W2f , 


(64) 


which can readily be checked numerically. Further, we also check numerically that 
introducing particle-hole excitations results in the following conformal weights with respect 
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to the ground state of the untwisted c = 2 theory 


A. 


A„ = - 
" 2 


1 — + 5mi^) — + - | 6mi6m2 


P 

2 2 
+ O + U!lU’2 + W2^) H-7 

3 ^ ' p — A 

+ + ^2 • 


{wi + W2y 


(65) 


It is also of interest to study what happens to these excited states when a twist (p is 
turned on. We therefore measure numerically the corresponding effective central charges, 
c« = 2-12 (A« + A„), as a function of the twist. Solving the corresponding BAE for sizes 
up to L ~ 100, we hnd the corresponding formula 


Ca{p) = 2-6 


p 


1-j [Sm^ + Smy) — [ 1 + - J 6mi6m2 


P 


- 8 [w^ + W 1 W 2 + W2^) - {wi + 102^ 


P 


6p (f 


^—h 2(ryi + 1 ^ 2 )^ — 12 (^Sf + 5]^ + + (^2 ) • (66) 


p — 4 TT 

An example is shown in £gure[^ where we measure numerically the central charge associated 
with the ground state (72) as a function of the twist, and which shows an excellent agreement 
with the formula (66). 

The exponents Aq,, Aq, are measured with respect to the untwisted model with central 
charge c = 2. Anticipating on the results of section 5.1.1 , the exponents of the Potts model 
should instead be measured with the twisted central charge (73). We shall therefore consider 
the shifted exponents ha, ha, 


Or, more precisely. 


ha T ha 2 


ha + ha — Aq, + Aq — 


p{p - 4) 


1 -^ + 5 mi^) — ^1 H—^ 6mi6m2 

+ X + W 1 W 2 + W2^) H- (wi + W 2 )" 

3 p — 4 

+ + 82 

■ ^ "(^^(—))-y^4, 


(67) 


( 68 ) 


2(p — 4) TT 

To disentangle ha from ha, we can further use the translational invariance of the system, 
implying that the TM commutes with the momentum operator. The latter can be written 
as a sum over Bethe integers. 


mi 


m 2 


p = 


E + E 2A 


. i=l 


k=l 


2tt 

T 


(69) 
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So we readily have, for the states of interest: 

21 
~3 


P 

“ L 


^ - Smij '^1 + f ^ “ *^"^2 ) ta 2 + - <^1 - <^2 


= 271 ^ ('^l' + “ <^1 ~ ^2 ~ - 61712102 ) 


(70) 


In the continuum limit, the ^^-dependent part is related to the conformal spin ha-ha, as 


ha-ha = 5 ^ + 5 ^ — ( 5 ;^ -62 — 6miWi — 61712102 ■ 


(71) 


5. Critical content of the colonring problem 

We will now restrict to the specihc states observed in the Potts spectrum. This will lead 
to analytic expressions for the central charge, the thermal and magnetic exponents, and the 
watermelon exponents describing the propagation of several FK clusters. Combining this 
information will hnally lead to the identihcation of the CFT describing vertex colourings of 
the triangular lattice. 


5.1. Ground state (k = 0) sector 

5.1.1. Central charge. As already explained, the ground state of the Potts model (namely, 
in the /c = 0 sector), can be found as an eigenstate in the zero magnetisation sector of the 
spin-one vertex model, at twist ip = 9. In terms of Bethe roots, it simply corresponds to 


{6mi,6m2) = (0,0), {wi,W2) = (0,0), tp = 9. (72) 


From equation ( |66[ ), we therefore directly read the corresponding central charge, 

24 

c = 2 


p(p - 4) 


(73) 


We have made an extensive numerical test of this formula by 

• Direct diagonalisation of the TM for small values of L = 6,9,12. In this case we were 
able to cover the whole interval Q G (2,4) in steps of size 10“^. 

• Numerically solving the BAE for sizes L = 6,9,12, 24,48, 96 with the right twist ip = 9. 
Except for L = 6, 9, we could solve these equations in the smaller interval Q G [2.469,4), 
again in steps of size 10“^. Note that for Q G (2, 55 ), the root conhgurations contain 
one or more 2-strings. Furthermore, the higher the value of L, the larger the number of 
2-strings that might appear in the conhguration. This technical complication prevented 
us from studying the interval (2, 2.469) for L > 12 by solving the BAE. Indeed, for 
L = 6,9,12 the results obtained from this method agree perfectly well with those 
obtained from the TM diagonalisation. 
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Potts model in Regime IV model in Regime IV 



Figure 7: (a) Values of the free energies fi as a function of Q for several values of the width: 
L = Q (red), L = 12 (blue), L = 24 (pink), L = 48 (orange), and L = 96 (violet). The 
dark-gray curve diverging at Q ~ 2.224 corresponds to L = 9. The black curve corresponds 
to Baxter’s exact result ( [l6| ). (b) Absolute difference |/buik — /Baxter| for the bulk free energy 
obtained from the three-parameter £t friQ) = /buik(Q) + ^c{Q) -|- ^ with = 6 (red), 
Lmin = 12 (blue), and Lmin = 24 (pink). 


In particular, for each value of Q G [2.469,4), we computed the free energy fiiQ) for 
L = 6,12,24,48,96. (fgiQ) was computed mainly for testing our methods.) We then 
performed the standard three-parameter £t based [cf. (35)] on the CFT Ansatz fiiQ) = 
/buik(Q) + l^c(Q) -|- using three consecutive values of the width, L = Lmin, 2Lmir,, 
to obtain the unknowns /buik(Q), c{Q) and A. As a precaution against hnite-size effects, we 
repeated the fits for three different values of Amin. 

The values obtained for /bulk agree very well with those given by Baxter’s exact 
expression for regime IV (16), even though we are also computing the former values for 
Q much less than Qq [cf. (19)], where the ground state of regime I is dominant. This vividly 
illustrates the fact—discussed at the end of section |2.4|— that regimes I and IV coexist, in 


the sense that the scaling levels of either of them are well-defined by analytic continuation 
on the other side of Qg. 

In hgure [^a) we depict the raw data for /l(Q), as well as Baxter’s free energy (16), 
and in panel (b) we show the discrepancies between /bulk and /Baxter (16) in semi-logarithmic 
scale. The agreement is very good. This result shows that the ground state of regime IV 
can be analytically continued deep inside the regime I. We believe that this agreement can 
be extended to the range 2 < Q < 4, and we have made extensive tests to verify this claim. 
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Indeed we have followed the free-energy level fi which is dominant at Q = 4 (i.e., the ground 
state of regime IV) by analytic continuation to smaller Q —using both direct diagonalisation 
of the TM for small L, and numerical BA computations for larger L —and found it to have 
the following features: 


• If L G 6N, the level /l follows Baxter’s g 2 closely down to some value Qmm{L) > 2 
where fiiQ) has a minimum. For Q < Qmin{.L)-i fi breaks away from g 2 . However, 
for L = 9 (see hgurej^a)), although the behaviour is similar, we hnd a divergence at 
Qdiv(9) ~ 2.244, rather than a minimum. Therefore we expect parity effects, so in our 
main hts we will only include even values of L. 

• There is good evidence that Qmin(-h) —)■ 2 as L —)■ cx) for L G 6N. We also conjecture 
that Qdiv(-h) —)• 2 as L —)■ cx) for L being an odd multiple of 3. 

• Right at Q = 2, we hnd that fi takes a hnite value that, at least for L = 6,9,12 
decreases as L increases. If we ht the corresponding /l( 2) values to the standard CFT 
3-parameter Ansatz, we obtain /buik(2) ~ —0.2570, which is still larger than Baxter’s 
value ^ 2 ( 2 ) ^ -0.4315231. 


These observations indicate that /^ —)■ (72 as L —)■ 00 on the interval Q G (2,4], but the 
convergence is not uniform. This feature is responsible of the ‘largest’ discrepancies observed 
for Q ^ R 5 . Indeed, for Q —2+ the central charge (73) diverges, signaling the break-down 
of regime IV. On the other hand, we have found no evidence of regime IV levels whose 
analytic continuation converge to g 2 for 0 < Q < 2 . Since g 2 is analytic on the whole interval 
0 < Q < 4, it is possible that it serves as the ground state of yet another regime (and a 
different CFT) for Q G [0,2), but further examination of this issue is beyond the scope of 
the present paper. 

A further noteworthy feature concerns the value of the analytic continuation of 
regime IV’s ground state to Q = 3. We hnd that /l( 3) = 0 exactly for any L. This 
should be compared with the simple fact that a triangular-lattice strip with cylindrical 
boundary conditions and width L a multiple of 3 is uniquely 3-colourable modulo global 
colour permutations (i.e. Zg{^, —1) = 3!). Therefore regime IV captures an essential feature 
of the colouring problem at Q = 3. 

Finally, in hgure [^b) we observe that the convergence of to g 2 is noticeable slower 
close to Q = 4. This is probably due to the usual logarithmic corrections at Q = 4, well- 
known in the Potts-model context. 

The reader might worry if all of this does not contradict the simple observation that 
Z = 0 for Q = 2, since a triangular lattice manifestly does not admit a proper vertex 2- 
colouring. As a matter of fact, there is no problem, because at Q = 2 massive cancellations 
of TM eigenvalues will occur j|] due to three distinct mechanisms: 

(i) Some of the eigenvalues are simply zero; 


f The same happens for other special values of Q, which are the Beraha numbers for cylindrical boundary 
conditions 


10 


and Q S N for toroidal boundary conditions 


12 
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Potts model in Regime IV model in Regime IV 



Figure 8: (a) Values of the central charge c obtained by performing the 2-parameter £t 
hiQ) - /Baxter (<5) = ^c{Q) + ^ for Lmin = 6 (red), Lmin = 12 (blue), Lmin = 24 (pink), 
and Ljnin = 48 (orange). The black curve shows the conjectured value (73). Note that the 
hts cover the interval (2,4) for Lmin = 6, while they cover the smaller interval [2.469,4) 
for the other values, (b) Relative difference between the estimated central charge and the 
conjecture (73) in semi-logarithmic scale. 


(ii) Other eigenvalues have a zero amplitude in the decomposition of the Markov trace 35 ; 

(iii) Some eigenvalues are equal up to a sign, and their combined amplitude summed over 
the various sectors of the TM is zero. 

The end result is that Z = 0 indeed, despite the existence of finite eigenvalues, and we have 
checked this explicitly for small L. 

Now that the agreement of /bulk with /Baxter has been hrmly established, we can obtain 
more precise estimates for the central charge c by performing hts with only two parameters, 
c{Q) and A, using the CFT Ansatz /l(Q) — /Baxter(Q) = ^ and the exact result 

(16). Accordingly, we htted our data for two consecutive values of L, namely L = Lmin, 2Lmin, 
for different values of Tmin- The results are depicted in figure [8j The agreement between 


the numerical estimates for the central charge and conjecture (73) is excellent, including for 


Q < Qo- There is indeed good evidence that c converges to (73) for all Q G (2,4], but once 
again the convergence is not uniform; note in particular that c —)■ — oo as Q —>■ 2+, whereas 
/Baxter stays Fnite at Q = 2. 

At Q = 3, we hnd that c = 0 exactly for any L, as expected from the total absence 
of hnite-size effects noticed above. This value c = 0 agrees with the obvious non-critical 
nature of the 3-colouring problem. On the other hand, close to Q = 4, the observed larger 
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L/3 



Figure 9: Configurations of the Bethe roots associated with the ground state in regime IV, as 
a function of the twist parameter. For (p < tt, it is described by a Fermi sea of ^ real roots 
(black dots) and a Fermi sea of ^ roots with imaginary part | (white dots). As —)■ 7r“, the 
real parts of the rightmost roots of each Fermi sea go to +oo. As p crosses tt, these two roots 
form one 2-string (black squares). The same process occurs once again as (p —?• 27r“, resulting 
in a conhguration with two degenerate 2-strings at p = 2 n. Although on the picture these 
are distinct, it should be kept in mind that the corresponding roots actually coincide two by 
two. 


discrepancies are again due to logarithmic corrections. 


5.1.2. Thermal exponent. It is now interesting to further consider the excited states in the 
k = 0 sector, and in particular the hrst of these, which in the continuum limit corresponds 
to the thermal field. Although the corresponding root conhguration is highly unstable under 
variations of the twist, we can nevertheless describe its analytic continuation at zero twist, 
which is made of 

• ^ — 2 real roots. 

• ^ — 2 roots with imaginary part |. 

• Two degenerate 2-strings, with zero real part and the same imaginary part. 


Let us now go back to the ground state conhguration (72), and try to follow it 
continuously as the twist is increased. The evolution of the Bethe roots from p = 0 
to p = 271 is depicted schematically in hgure It is seen that the root conhgurations 
encounters qualitative changes, namely, the formation of 2-strings, as p crosses multiples of 
71. More precisely, the root conhguration associated with the thermal exponent at p = t) 
actually corresponds to the continuation of the ground state conhguration to p = 27i. The 
corresponding conformal weight can therefore be obtained straightforwardly by analytic 
continuation of the conformal weight for the twisted ground state (this is after all exactly 
the same as what was observed in regime I, see section 4.2) [§] 


§ More generally, we expect that all root configurations with 2-strings can be understood as the continuation 
of ‘standard’ root configurati ons to larger values of the twist parameter. This explains a posteriori our 
comment at the end of section [TSl and we shall not comment anymore on the presence of such configurations 
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Potts model in Regime IV (p=0 



Q 



Potts model in Regime IV <p=0 



Q 

(b) 


Figure 10: Fits for the subdominant free energy in regime IV at zero twist, (a) Values of the 
effective central charge Cefr obtained by performing the 2 -parameter £t / 2 ,l(Q) — /Baxter(< 5 ) = 
I^Ceff(Q) -h for Lmin = 12 (red), L^m = 24 (blue), = 48 (pink), L^m = 96 (orange), 
and L = 192 (violet). The black curve shows the conjectured value (74). (b) Relative 
difference between the estimated central charge and the conjecture (74) in semi-logarithmic 
scale. 


To check this assumption, we have solved numerically the zero-twist BAE associated 
with the first excited state continued to = 0 (remember that the true level present in 
the Potts model, namely at a twist 99 = 6 *, corresponds to a highly unstable configuration 
of roots, whose numerical study does not appear feasible) for sizes up to L = 384 in the 
interval Q G For each value of Q, we computed the corresponding free energy 

/ 2 ,l(Q) for L = 12,24,48,96,192, and 384. We then performed the standard three- 
parameter £t f 2 ,L{Q) = /buik(< 5 ) + f^CesiQ) + for three consecutive values of the width 
L = Lmin, 2Tmin, 4Lmin- As for the leading eigenvalue, the values obtained for /bulk agree very 
well with those given by Baxter’s exact expression for regime IV (16). A better estimate for 
Ceff is obtained by using the 2-parameter CFT Ansatz / 2 ,l(Q) - /Baxter(Q) = i^Ces{Q) + -A, 


and two consecutive values of L = Luiin, 2 Lmin. 


The results are depicted in hgure 10, and 


agree with the conjecture 


Ceff — —22 — 


96 

p — 4 ’ 


(74) 


that is, precisely the expression obtained from ( | 66 [ ) at Smi = 61712 = wi = W2 = 6f = 6f = 0 
for a twist = 27r. The agreement between the numerical estimates for the effective central 


in the other Potts sectors. 
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charge and conjecture (74) is excellent, especially for values of Q far away from 4. Close 
to Q = 4, we should expect logarithmic corrections as those shown in hgure implying a 
slower convergence towards (74). 

From there we assume, as in regime I, that the subleading eigenlevels in the k = 0 sector 
are obtained from the ground state level by taking twists of the form (p = 0 + 2 nj, leading 
to effective central charges 

24 9 

Ceffj = 2 - — -(1 + jp) , 


P{P - 4) 

and therefore to conformal weights of the form 

1 2 j 

hj + hj = (c - Ceff,i) = (2 + jp) . 

Among these, the two least irrelevant correspond to j = 1 and j = —1, namely 

4 , , ^ 12 


X'j’ — h_i -|- h_\ — 2 -|- 


p — 4 


Xrp — h\ hi — 2 -|- 


p — 4 


(75) 


(76) 


(77) 


5.1.3. Magnetic exponent. As explained in section 3.4, the magnetic exponent is associated 
to the scaling of the leading TM eigenvalue in sector k = 0, with twist f “ f • From our 
observations in section 5.4[ we know that this eigenvalue can be associated with the quantum 
numbers 

((5m,,im2) = (0,0), (wi.wz) = • ‘7 = 6. 

(79) 

(80) 


2 ’ 2 / 

and therefore with an effective central charge 

Ceff 0 . 

From there, we directly read off the magnetic exponent, 

1 7 ^ 1 2 


5.2. Even k sectors. 

For instance we study the sector k = 2, for Q 
are of the form fi=p log Aj, where A j are the 


/o 

= 0.356989 

{ 6 mi, 61712 ) 

/l=/2 

= 0.343760 

{ 6 mi, 61712 ) 

/s = 

= 0.321961 

(5mi, 61712 ) 

/s = /e 

= 0.316346 

{ 6 mi, 61712 ) 

h = fs 

= 0.314850 

{ 6 mi, 61712 ) 

/g = /lo 

= 0.309532 

{ 6 mi, 61712 ) 

/ii = /12 

= 0.304716 

(5mi, 61712 ) 

/l3 = /l4 

= 0.302223 

{ 6 mi, 61712 ) 


= 3.99, and L = 9. The eigenenergies given 
TM eigenvalues: 

(1.1) ,(w;i,'«;2) = (0,0),(p = 0 

(1.1) , (wi, W 2 ) = (±1, 0), p = =F7r 

( 1 . 1 ) , (wi, W 2 ) = (± 1 , 0 ), <p = 0 

(2, 2) + one 2-string 

(1.1) , (wi, wa) = (0, 0), = 1, p = 0 

(1.1) , {wi, W 2 ) = (±1, ±1), (p = 0 

(1,1), (wi, W 2 ) = (±2, 0), p = =F27r 

(1,1), (wi, W 2 ) = (±1, 0), = 1, p = 0 
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(when no mention is made of the (5i, 62 numbers, these are equal to zero). 

From this example as well as further examination of the k = 4,6,8 cases, we see that 
the least excited states making up the even k >2 sectors are associated with 


6 mi = 6 m 2 = -, wi,W 2 GZ, ip = nn, neZ, 
From we hnd the corresponding conformal spins 


(81) 


h-h = --k {wi + W 2 ) , 


(82) 


which indeed assume integer values, while the conformal weights are found to be 


h 

h 


[k{p — 4) — {2n + 2{wi + W 2 ))pf — 16 1 


+ — {3kn + {w 2 - wi) ) 


16p{p - 4) '12 

[k{p — 4) — {2n + 2{wi + W 2 ))p]^ — 16 
16p{p — 4) 

+ ^ {3k{n + 2wi + 2 W 2 ) + {w 2 - wif) . 


(83) 


(84) 


5.3. Odd k > 3 sectors. 


We here specialize to odd sectors with k > 3, as the k = 1 case is a bit particular, as we 
shall see in the next section. Similarly, for k = 3, Q = 3.99, and L = 9, the lowest-lying 
levels are the following: 


/i = /2 = 0.316285 

/3 = /4 = 0.313692 

/5 = /6 = 0.290629 

/7 = /8 = 0.288732 


{Smi, 6m2) = (2,1), (tci, ^ 2 ) 
{6mi, 6m2) = (2,1), (tci, ^ 2 ) 

{ 6 mi, 61712 ) = ( 2 , 1 ), (wi, w; 2 ) 
{ 6 mi, 61712 ) = (2,1), {wi,W 2 ) 


= (d=i, 0 ),(p = 0 
= (±|,0),(p = Tf 
= (±|,0 ),(p = 0,5^ = l 
= (±i,=Fl) ,9? = 


For k = 5, Q = 3.99, and L 

h = f 2 = 0.239272 

/s = /4 = 0.229800 

h = h = 0.224099 

h = U = 0.222247, 

/g = /lo = 0.220503 


9, we hnd 

{6mi, 61712) = ( 3 , 2 ), (wi, W2) 
{6mi, 6702) = ( 3 , 2 ), (tci, ^2) 
{6mi, 61712) = ( 3 , 2 ), (wi, W2) 
{6mi, 61712) = ( 3 , 2 ), {wi, W2) 
{6mi, 61712) = ( 3 , 2 ), {wi, W2) 


{0,±\),p = Tf 
(0,±|),(p = Tf 

(Tl,±|),(p = 

(o,±i),^ = Tf .y = 


In the latter case states corresponding to (p = 0, for instance { 6 mi, 61712 ) = (3, 2), (tci, W 2 ) = 
(0, ±|), 93 = 0 are also found, but lie deeper in the Potts spectrum. 

From these examples, as well as further examination of the k = 7 and k = 9 cases, we 
see that the least excited states making up the odd k > 1 sectors are associated with 


k -\- 7 
6 m, = — 


6m2 


k -1 

2 


k - 1 

Wi = Vi-\ --— 


W2 = V2 + 


k + l 
4 


V? 


2 mi 


( 85 ) 
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where Vi, V 2 and n take integer valnes. 

From we find the corresponding conformal spins 

h-h=|(l-/c^ + 2{v 2 - vi) - 2k{vi + ^ 2 )) , 

which indeed assnmes integer values, while the conformal weights are found to be 

^ _ [Hp - 4) - (/c + 4f + 2{vi + V2))p] ^ - 16 

16p{p — 4) 

+ ^ ( 6 n + (^2 - + 2 )^) 

- _ [k{p - 4) - (A; + 4f + 2{vi + V2))p] ^ - 16 
16p{p — 4) 

+ ^ + 6(n + kvi + kv 2 ) + (^2 - Ui - 1)^) 


( 86 ) 


(87) 


( 88 ) 


5 . 4 . k = 1 sectors. 


As we mentioned in section 3.4, the k = 1 case is a bit particular, in that there are two 
ways to have one propagating cluster in the vertical direction, namely it may or may not be 
allowed to wrap around the cylinder. Let us study these two cases independently: 


5 . 4 . 1 . Cluster allowed to wrap around the cylinder. Forbidding loops on the kagome lattice 
to wrap around the cylinder amounts to giving such loops a weight 0, and therefore to fixing 
a twist 2 = f • 

Indeed, the largest eigenvalue of the FK TM in this sector is found in the zero- 
magnetisation sector of the vertex model at twist p = tt. Alternatively, it can be obtained 
at zero twist from the following set of Bethe parameters 


(imj.Smz) = ( 0 , 0 ), (wj.wz) 



(89) 


5 . 4 .S. Cluster not allowed to wrap around the cylinder. Forbidding the propagating cluster 
to wrap around the cylinder amounts to fixing one through-line in the 0{n) model. Therefore 
we expect to find the corresponding leading eigenvalues in the sector of magnetisation m = 1. 

This is the case as, for instance, we find for Q = 3.99 and L = 6 the following set of 
leading eigenenergies: The first dominant states form complex quartet 

/i = 0.302084 + 0.3553561 
/2 = 0.302084 -0.3553561 
fs = 0.302084 + 0.3427751 
U = 0.302084 -0.3427751 





40 


Q-colourings of the triangular lattice 


and the next states are given by 

h = h = 0.272659 
h = h = 0.269606 
/g = 0.227715 


(5mi, 8 m2) = (2, -1), (^1,^2) = (±|, O) , 9? = 0 
(5mi, 5 m2) = (2,1)+ one 2-string, = 0 

( 5 mi, 5 m 2 ) = (2, -1), (^1,^2) = (0, 0), 92 = 0, 

5 + = 57 = i 


The complex qnartet, as we nnderstand it, is associated to very pecnliar root 
conhgnrations. From observations at L = 6 and L = 9, we see that the corresponding 
root conhgnrations are made of: 

• ^ — 1 roots with imaginary parts very close to 0 . 

• ^ — 1 roots with imaginary part very close to |. 

• One root with imaginary part very close to |. 

Snch conhgnrations are rather odd, and ont of the scope of onr analytical solntion. In 
addition, we were nnfortnnately nnable to solve the corresponding BAE nnmerically for 
large system sizes, and therefore we could not access the corresponding conformal weights. 


5.5. Identification of the CFT 


The central charge (73) is the analytic continuation to any real p > 4 of that of the 
S' 3 -invariant minimal models studied by Fateev and Zamolodchikov 63, 64 , which can 
alternatively be formulated as the conformal cosets 

Sf/(2)4 X Sf/(2)p_e 

dehned for integer values of p > 6 . 

The corresponding sets of conformal weights are known 6^ 64 to be parametrised by 
two integers, l<r<p — 5, l<s<p — 2, and read 

[rp — s(p — 4)]^ 


hr.fi — 


16 ^ 1 — cos*^ ('^^) 


16p(p - 4) ' 12 

In addition there is another series of primary helds, with conformal weights 

1 -|- sin^ 

H- :: -1" + 2(5r,2^s,l • 


hrg hj. g 


(91) 


(92) 


Both (91) and (92) can be written under the form hr^si.h'j.g) = where 

Shr^s encodes the sector of the associated S '3 symmetry [^, and can be understood in the 
framework of a generalized CG construction as a boundary twist term (this is analogous 
to the Neveu-Schwartz and Ramond sectors in supersymmetric CFTs, encoding the periodic 
or antiperiodic boundary conditions for the fermionic field |^). 

Let us now compare the conformal spectrum identified in regime IV with the exponents 

(P/dH- 
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5.5.1. Thermal exponents. In section 5.1.2, we have identified a series of subleading 


exponents, hj + hj, given by (76). These exponents have zero momentum (in agreement with 
the fact that thermal fields should not break the translational symmetry), and therefore 


hj = hj = —- 

p-A 


(2 + jp) . 


(93) 


These are readily identified as the dimensions in (91), with r = 4j + 1, and s = 1. 


5.5.2. Magnetic exponent. In section 5.1.3, we have identified the magnetic exponent as 


(80). The latter cannot be written in terms of the table {h'rs) with fixed, p-independent 
values of r, s. Instead, we find that it coincides with xh = h^^s + with 


p — A p 

r = -, s = - . 

2 ’ 2 


(94) 


The fact that these indices are p-dependent and half-integer is not too surprising. Indeed, 
the magnetic exponent for the usual critical Q-state Potts model has similar features, since 
it can be written h p-i p+i , where p here denotes the minimal model indexj 


5.5.3. Even k sectors. In section 5.2, we have identified the leading conformal weights in 


the sectors with even k > 2 as given by eqs. (83)/ ([84)). Note that the expressions for h and 


h can be recovered from each other by se 


ling 


Wl,W2) = 


-Wi, 


-W 2 ), n' = n + 2wi + 2 w 2 , 
so we only focus on the former, namely (83). The first term takes the form 


with r = 2n + 2{wi + W 2 ), and s = k. Focusing on the most relevant state, namely that with 
n = Wi = W 2 = t), the second term vanishes, and we recover for fc = 0 (mod 4) an exact 


correspondence with the exponents (91), for which the term 5hr^s is zero. This state has zero 
conformal spin and scaling dimension 


Xk = h + h = 


(p — 4) 

8 p 


pip - 4 )' 


(95) 


5 . 5 . 4 . Odd k sectors. In section |5.3[ we have identified the leading conformal weights 
in the sectors with even k > 2 as given by eqs. (87)/(88). As in the k even case, the 


expressions for h and h can be recovered from each other by a relabeling of the indices, namely 




-wi. 


'W 2 ),n' = n+k{wi+W 2 ), so once again we only focus on the former, namely 
(87). The first term takes the form ? with r = k+^+2{vi+V2) = ^+2{wi+W2), 

and s = k. In particular, we notice that fractional values of r are allowed in this case. It 
is also noticeable that contrarily to the case of k even, no well-definite set of values of the 
parameters Ui,U 2 ,u seems to give rive to a least irrelevant state. 


This result can however also be written as hi g, with p-independent indices. 


























Q-colourings of the triangular lattice 
5.6. Physical conclusions 


42 


The Q-colouring problem of the triangular lattice G is dehned sensu stricto only for Q G N. 
We have studied here its analytic continuation, in the form of the chromatic polynomial 
Xg{Q), defined for Q G M. Within the critical range 0 < Q < 4, two distinct regimes (that 


we have called regimes I and IV) were identified 17 , and we have shown that the ground 
state and scaling levels of either regime are well-dehned regardless of considerations 18 
about which ground state is dominant in absolute terms. In particular, regime I can be 
dehned for 0 < Q < 4, and regime IV makes sense for (at least) 2 < Q < 4. 


We now substantiate our claims—hrst made in section 2.4) — that regime IV is in fact 
the correct analytic continuation of the colouring problem. To this end, we show that this 
regime has the correct physical properties of a colouring problem for Q = 4 and Q = 3, 
and that its break-down in the limit Q —)■ 2+ can again be motivated by properties of the 
colouring problem. 

The Q = 4 vertex colourings of the triangular lattice are equivalent!^ to vertex 3- 


colourings of the kagome lattice and to edge 3-colourings of the honeycomb lattice 68,69 


The latter problem has been studied in its own right pi8j and, more generally, as the special 


case n = 2 of the 0(n)-type fully-packed loop (FPL) model on the hexagonal lattice 70,71 
The ground state free energy of these problems agree [T7l|^[70] and is given by (16), and 


since the formulation in terms of 3-edge colourings of the honeycomb lattice involves only 
real, non-negative Boltzmann weights [^, the Perron-Frobenius theorem ensures that (16) 
correctly gives the entropy of vertex 4-colourings of the triangular lattice. 


A two-dimensional height mapping for the 0(?7,) FPL model was given in 71 , and it was 


shown that the continuum limit is described by a Liouville field theory (or two-component 
CG) involving two compactified bosons with an n-dependent radius. For n = 2 this reduces 


to two free bosons, so the central charge is c = 2 in agreement with our result (73) for Q = 4 
{p —)■ cxd). a slightly different two-dimensional height mapping was given directly for the 
4-colouring problem in 69 . The critical exponent rj = 2 xh = | was also derived and checked 


numerically in 69 , in agreement with our result (80) that gives xh = h for Q = 4. This 


provides a couple of convincing checks that our continuum description is indeed that of the 
combinatorial colouring problem]^ 

we claim 


We next turn to the case Q = 3. Even though this is less than Qq [see (19) 


that regime IV still correctly describes the physical colouring problem. It is easy to see that 
each of the three sublattices must have a constant colour, so Z = 3! for any size L G 3N, 


and this is compatible with our observation, made in section 4.4, that /l(3) = 0 . A further 


^ Two models are here said to be equivalent if their set of states are in bijection—up to an appropriate 
choice of boundary conditions, such as choosing an appropriate (e.g. even) size of the lattice in any periodic 
direction and/or fixing the colour of a degree of freedom in one of the models—and if the corresponding 
Boltzmann weights are equal, up to any overall factor. Also, the equivalences mentioned require the lattices 
to be planar, and the situation on the torus might involve further subtleties (see 67 and references therein]). 


Note however the watermelons exponents of 71 do not compare directly to ours, since these observables 


are different from our spin-one watermelon operators. 
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confirmation comes from the magnetic exponent. Since there is a perfect ordering on each 
sublattice, the (staggered) two-point function is constant, so that xh = 0. This indeed agrees 


with (73) for Q = 3 (p = 6). 


Finally, consider the limit Q —?• 2+. It is obviously not possible to colour an elementary 
triangular face with only two colours, but this does not imply that the ground-state 
free energy must diverge in this limit, because the cancellation mechanisms mentioned in 


section 4.4 are at play at Q = 2. Indeed the regime IV expression (16) stays hnite in this limit. 


However, the central charge (73), and all critical exponents derived above contain a pole at 


Q = 2 (p = 4), which is concomitant with the break-down of the colouring interpretation. 
On the other hand, it is possible to shift the ground-state energy by an inhnite amount, 
so as to obtain a well-dehned Q = 2 theory in which each elementary face of the lattice 
contains precisely one frustrated edge (i.e., that is adjacent to two equally coloured spins). 
This dehnes the so-called zero-temperature Ising antiferromagnet, which possesses a critical 
ground-state ensemble corresponding to a c = 1 CFT of one free boson 72 -74 . It seems an 


appealing idea that this CFT should be (inhnitely deeply) buried inside regime IV, but we 
have not investigated this question in details. 

Further support that the critical colouring problem is described by regime IV for 
2 < Q < 4 is provided by the RSOS representation to be described in section 

It is worth noticing that the results found in this paper favour Baxter’s value [T^ for Qq 


[cf. (19)] against the alternative value 2 -|- suggested in [1^ section 6]. Furthermore, this 
conclusion may be a signal that Baxter’s eigenvalue gi ( p)5| plays no role inside the critical 
region of the phase diagram in the complex Q-plane of the triangular-lattice chromatic 
polynomial (as it was investigated in |^). Notice that even though gi is not dehned for any 
real Q G (0,4), it is well dehned when Im (Q) ^ 0. 

Before closing this section, we notice that our results for regime IV also have implications 
for the phase diagram in the (Q,n) plane, outside the chromatic line v = —1. According 


to (77), the thermal perturbation of the ground state in regime IV is irrelevant for all 
Q < 4. Therefore the part of the chromatic line where regime IV is dominant—that is, 
Qo < Q < 4—is attractive under the renormalisation group. We conclude that the small 
‘bubble’ in hgure which is delimited by the bifurcation point T, the two branches of 
and the vertical line Q = 4, is in the universality class of regime IV. 


6. Connection with RSOS models and anyonic chains 


In section 2^ we discussed that at the Beraha numbers, Q = 3^, all of the unphysical regime I 
cancels out of the partition function. What remains is a purely local model encompassing 
regime IV. This local model is conveniently formulated as a restricted solid-on-solid (RSOS) 
height model, which has independent interest because of its relation with anyonic spin chains. 













Q-colourings of the triangular lattice 

6.1. The RSOS model associated with the Izergin-Korepin vertex model 


44 


In section we have explained how the colouring problem on the triangular lattice could 
be mapped onto an appropriately twisted version of the IK model, which can be viewed 
as a spin -1 vertex model; the degrees of freedom on each lattice edge are vectors in spin -1 
representations of the quantum group t/q(s/ 2 ), where q = e *2 = e*p. These representations 
can at hrst stage be understood in the same way as classical SU{2) representations, and the 
interaction between two edges at a vertex of the square lattice can therefore be decomposed 
into projectors on the spin 0 , spin 1 and spin 2 representations, namely 



oc 



sin {9 + I) 

sin (e - f) 



cos (f + I) 


(a) 


cd 

ab 


(96) 


(for convenience the lattice has been rotated by 45° with respect to the figures in section]^. 

As we will now explain, this allows for a reintepretation of the model in terms of integer 
heights living on the faces of the square lattice (viz. vertices of the dual square lattice, shown 
in green in the following hgure), namely in terms of an interaction-round-a-face (IRF) or 
solid-on-solid (SOS) model. The construction is adapted from the seminal work of Pasquier 
in the case of spin-| models 
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, and can be summed up as follows (for more details see 


58 



the ji and j[ living on the nodes of the square lattice are (half)-integers encoding the 
representations of 17q(sZ2), and which must thus be compatible with the f/q(s/ 2 ) fusion rules: 
for instance, if representations ji and jj+i are separated by a lattice edge (carrying a spin -1 
representation), jj+i must belong to the fusion product j* x 1 . 

The action on the projectors Pq, Pi, P 2 on a given face of the dual lattice, specified by 
a set of heights ji_i, jj, jj+i, j', is worked out by a change of basis which makes use of the 
quantum analog of Wigner’s 6 j symbols and can be represented as follows : 


1 1 


ji+l 


E 


ji—l 1 ji 
1 ji+l J 


1 



ji 


1 


ji+l 


(97) 


3i—l 


Ji 


Ji-1 
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and similarly for the upper part of the face. Combining it all, the action of a projector Pj 
reads 



We now specify to the case of q equal to a root of unity, namely p integer, where the 
fusion rules of the quantum group truncate, leaving us with a restricted set of possible 

representations 


j = 




(99) 


The resulting IRF model then becomes an RSOS model, with a hnite-dimensional set of 
degrees of freedom per lattice site. The corresponding TM can therefore be computed, and 
we expect that its spectrum is related to that of the original vertex model. More precisely, 
it was observed in 58 that the leading eigenvalues of the RSOS model could be obtained 
from those of the vertex model in the m = 0 sector, for values of the twist p = 6,36 ,.... 
The ground state of the RSOS model is in particular obtained as the ground state in the 
m = 0 sector, for p = 6. Taking for the spectral parameter the special value = f ~ | 
ensuring the correspondence with the Potts model, this state is precisely the ground state 
of the Potts model in regime IV. 

Let us check this by comparing the spectrum of the two models (colouring model and 
RSOS) for hnite-size systems. 

We start with, for instance, p = 36, which corresponds to Q ~ 3.969616 > Qq. Among 
the eigenvalues of the RSOS TM for L = 6, we hnd that indeed: 


• The ground state corresponds to the ground state of the Potts model in the k = 0 sector. 

• The 11**^ and 12**^ eigenvalues correspond to the leading doublet of the Potts model in 
the sector k = 1, where the cluster is allowed to wrap around the torus. 

• The 42**^ eigenvalue corresponds to the leading eigenvalue in the sector with k = 2. 

• A very low-lying eigenvalue corresponds to the hrst excited state in the k = 0 sector of 
the Potts model. 

• Noticeably, no eigenvalue of the RSOS model was found to correspond to the leading 
quartet of the k = 1 Potts sector, where the cluster is not allowed to wrap around the 
torus. 

The conclusion we can draw from these observations is that, apart from the very peculiar 
non-wrapping k = 1 sector, the leading states of the Potts model in regime IV are all 
recovered by the RSOS model. Interestingly, the conclusion still holds when we take smaller 
values of p, for instance p = 6 or p = 7, for which Q < Qq: while the FK representation of the 
Potts model, and equivalently, the twisted IK vertex model are then dominated by regime I, 
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the ground state of the RSOS model is still that associated to regime IV (in particular, 
the corresponding free energy is zero when Q = 3). In other terms, the RSOS model, as 
the colouring problem, is described by regime IV for any value of 3 < Q < 4. While it 
seems legitimate, regarding our previous observations for 2 < Q < 3, to expect that the 
correspondence still holds down to Q = 2+, comparison of the Potts and RSOS eigenspectra 
for p = 5 (Q — 2.61803) shows that this might actually not be the case. Indeed, what we 
identihed as the continuation of the ground-state energy of regime IV for this value of Q was 
not found in the spectrum of the RSOS model, neither for L = 6, nor for L = 12. 

Therefore, we can conclude that the RSOS model is, for p > 6 and for a certain range of 
values of the spectral parameter, described by the OFT governing the physics in regime IV, 
namely the coset 

As we shall see in the next section, these conclusions can actually be related to recent 
observations relative to the physics of spin-1 anyonic chains (20). 


6.2. Spin-1 anyonic chains, brief introduction 


The RSOS model considered in the previous section can be viewed as the 1-|-1-dimensional 
formulation of a spin-1 anyonic chain, which is a deformation of the usual Heisenberg chain for 
SU{2) spins, well suited for the description of particles with exotic topological properties, the 


so-called non-Abelian anyons (see [^, or 75 for a more general and thorough introduction). 

For a given integer k, SU{2),^ anyons can be thought of as angular momenta indexing 
the representations of f/q(s/ 2 ), with q = eV = e*^, and therefore restricted to 


j = 0,t 1,... 


( 101 ) 


A spin-1 chain of SU{2)^ anyons can be represented as the following ‘fusion tree’ 


111111 


ji 32 3z 


where as usual the dashed lines indicate for periodic boundary conditions, and as in the 
previous section the anyonic indices ji must be compatible with the SU{f2),^ (f/q(sZ 2 )) fusion 
rules. As in spin chains with quadratic interactions, the Hamiltonian can be written as 

H = (102) 

i 

where encodes the projection of the product of two j = 1 representations over the 

decomposition 


1x1 = 0 + 1 + 2. 


( 103 ) 
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= cos92,i (-P2)i,i+1 


(104) 


The connection with the fusion tree representation above is worked out by changes of basis 


which rearrange the order in which the anyons are fused, which we already described in (97). 


The fusion tree also provides an easy means of computing the dimension of the RSOS 
TM. The adjacency matrix (in the target space of angular momenta) takes the form of an 
M X M matrix with rows and columns indexed by 2j + 1, shown here for M = k + 1 = 6 : 


/ 0 1 0 0 0 0 \ 
1110 0 0 


A = 


0 1 1 
0 0 1 
0 0 0 


1 0 0 
1 1 0 
111 


(105) 


\ 0 0 0 0 1 0 / 

Note the O’s in the upper left and lower right corner. The eigenvalues are 


Ak = 4cos^ 


kir 


- 1 


(106) 


with k = 1,2,...,M, so the dimension of the TM for a periodic chain of length L is 
At ■ This is dominated by the k = 1 term, so the asymptotic behaviour is 

{Q — 1)'^, the same as the initial colouring TM. 


6.3. Spin-1 anyonic chains: CFT results 

The phase diagram of spin-1 anyonic chains was recently investigated in 20 ,|]giving 

evidence for a wealth of gapped or critical phases as the angle 02,i is varied. An integrable TM 
and a one-dimensional quantum Hamiltonian are usually connected via the very anisotropic 
limit M —7- 0 of the spectral parameter. This allows one to relate the lowest Hamiltonian 
energy levels to the leading TM eigenvalues for extended regions of the spectral parameter. 
In the present case, things are slightly different, since the regime IV was observed to dominate 
the RSOS spectrum for a range of spectral parameters centered around the special value 
n = I — |. This range does not contain the anisotropic limit n —)■ 0, and in this sense, there 
is no straightforward relation to any properly defined Hamiltonian. It is however instructive 
to consider the Hamiltonian H derived from the m —>■ 0 limit of the integrable TM, which, 
in the vertex formulation, is dominated by either regimes I or HI (depending on the overall 
sign of H) for the values of 9 considered in this paper (namely 0 < 0 < |). As shall be 


reported in more detail in 76 , we observe that the eigenlevels corresponding to regime I are 


f Further details have been worked out by two of us and will be presented elsewhere 76 . 
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actually absent from the RSOS formulation, and that the associated RSOS Hamiltonian is 
instead dominated by those of regime IV. The associated value of the parameter 62^1 is 


02,1 = vr — arctan 


27r 47r ’ 

cos - cos 


(107) 


and lies precisely inside an extended phase identihed in 20 as being described by the 


conformal coset ( 100 ). 


7. Conclusion 

In this paper we have revisited the problem of Q-colourings of the triangular lattice hrst 
investigated by Baxter 17,18 . Using Bethe Ansatz techniques, this author found three 


different analytical expressions gi, g 2 , gs for the free energy in the thermodynamical limit. 
The latter two describe the critical range 0 < Q < 4, with g 2 being dominant for Q G (Q 05 4] 
and g 3 being dominant for Q G [0, Qq), where Qq = 3.819 671 ■ 


18 


We have gone beyond these results by determining the conformal held theories (CFTs) 
describing the scaling levels above (regime I) and above g 2 (regime IV). This was done 
by mapping the model to a spin-one vertex model on the square lattice, again amenable to 
Bethe Ansatz analysis. In particular we have found analytical expressions for the central 
charge and various critical exponents for both regimes. 


While regime I is well-known from the spin-one loop model on the square lattice 57 


regime IV is new and emerges because of the particular way the vertex model must be 
twisted in order to ensure the equivalence with the initial colouring problem. Accordingly 
we have dedicated the bulk of the paper to regime IV. It corresponds to a system of two 
compact bosons and we have identihed the corresponding coset theory. We have also related 


its RSOS representation to a recently studied spin-one anyonic chain 20 


One of our main conclusions is that regime IV possesses the correct physical properties 
to describe an actual colouring problem. By analytic continuation it can be extended to the 
range Q G (2,4], although the ground state g 2 is only dominant for Q G {Qo,4\. In this 
sense the two CFTs coexist throughout the critical range. We believe it to be crucial that 
regime IV describes features of the actual colouring problems at both Q = A and Q = 3, 
breaks down for Q —>■ 2 + as expected from the impossibility of two-colouring a triangular 
lattice, and that its restriction in terms of a local RSOS height model maintains the critical 
exponents of Q-colourings at the Beraha numbers (|^ for all integer p > Q. 

We have however left for future work the interesting question of what is the correct 
scaling theory above g 2 for Q G [0, 2). 

The case Q = 4 of regime IV is equivalent to several other models, including fully- 
packed loops (FPL) on the honeycomb lattice with fugacity n = 2. The continuum limit of 
both these models consists simply of two free bosons. However, while a Coulomb gas (CG) 


deformation of the FPL model to n < 2 is well-known 71 , the CG description of the Q < 4 


colouring problem is presently unknown. Identifying it presents another challenge for further 
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investigations, and might eventually permit one to rederive the critical exponents given here 
much easier (albeit less rigorously), without any use of the Bethe Ansatz 


71 


Finally, it would be interesting to study boundary extensions of the Q-colouring problem, 
by replacing the geometry of the cylinder (torus) with that of an 


along the lines of 42 


annulus. In particular, one could expect from jW,78 that restricting the number of colours 
allowed on the boundary would lead to sets of continuously varying critical exponents. 
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